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INTRODUCTION 


The  Thirteenth  DoD  Tri-Service  Review  Conference  on  Atmospheric  Transmission  Models 
was  held  at  the  Geophysics  Laboratory,  Hanscom  AFB,  Massachusetts  on  5-6  June  1990.  The 
purpose  of  the  meeting  was  to  review  progress  in  the  modeling  of  radiation  propagating  through 
the  earth’s  atmosphere,  identify  deficiencies  in  these  models,  and  make  recommendations  for 
improvements. 

Approximately  100  scientists  and  engineers,  representing  DoD,  other  government  agencies, 
industry,  and  the  academic  community  were  in  attendance.  The  agenda  consisted  of  thirty-one 
papers  with  sessions  on  turbulence,  atmospheric  propagation  models,  aerosols,  UV  modeling, 
and  LIDAR. 

This  proceedings  volume  summarizes  the  technical  presentations  at  the  conference.  The 
main  part  of  the  report  consists  of  the  abstracts  and  copies  of  the  viewgraphs  or  slides  and  other 
material  as  provided  by  the  authors  for  the  presentations.  The  Appendix  includes  the  original 
cril  for  papers  and  invitation  to  the  meeting,  a  copy  of  the  Agenda  for  the  meeting,  a  list  of 
the  attendees,  and  an  author  index. 

A  special  thanks  is  extended  to  Ronald  G.  Isaacs,  Betty  Stenhouse  and  Atmospheric  and 
Environmental  Research,  Inc. 

Francis  X.  Kneizys 
Leonard  W.  Abreu 


STATUS  OF  THE  LOWTRAN  AND  MODTRAN  MODELS 

L.W.  Abrcu,  F.X.  Kncizys,  G.P.  Anderson,  E.P.  Shcttlc*,  and  J.H.  Chetwynd 
Geophysics  Laboratory/OPE,  Hanscom  AFB,  MA  01731 

LOWTRAN  7  (1988,  AFGL-TR-88-0177)  calculates  atmospheric  transmittance  and  background 
radiance  for  any  given  atmospheric  path  using  spherical  refractive  geometry.  The  calculations  arc  done 
at  low  spectral  resolution  (20  cm  '  full  width-half  maximum).  The  effects  of  single  solar/lunar  scattering 
and  multiply  scattered  thermal  and  solar  radiation  are  included.  Enhancements  and  modifications  to 
LOWTRAN  will  be  discussed. 

MODTRAN  (presently  beta-test  version)  is  a  2  parameter  band  model  (P  &  T)  code  with  moderate 
spectral  resolution  (2  cm  *  full  width-half  maximum).  The  MODTRAN  band  model  parameters  were 
developed  by  utilizing  the  HITRAN  data  base.  Full  compatibility  with  LOWTRAN  7  is  maintained  and 
the  code  can  be  degraded  to  between  2  and  50  wavenumber  spectral  resolution. 

‘Now  at  Naval  Rc.scarch  Laboratory  (NRL),  Code  6520,  Washington,  D.C.  20375 
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STATUS  OF  THE  LOWTTUIM  AiND  MO[yrRi\N  MODELS 


L.W.  Abreu,  F.X.  Kneizys,  G.P.  Anderson, 
E.P.  She t tie  and  J.H.  Chetwynd 
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MODELING  SOLAR  AND  INFRARED  RADIATION  FIELDS  FOR  BTI/SWOE 


J.R.  Hummel 

SPARTA.  Inc.,  24  Hartwell  Avenue,  Lexington,  MA  02173 

The  Balanced  Technology  Initiative  (BTI)  on  Smart  Weapons  Operability  Enhancement  (SWOE)  has 
a  goal  to  model  the  radiant  field  from  complex  natural  backgrounds.  In  order  to  achieve  this  goal,  one 
must  be  able  to  model  the  thermal  structure  of  the  background.  Two  of  the  inputs  to  the  energy  budget 
of  the  background  arc  solar  and  infrared  radiation. 

A  review  of  approaches  to  calculate  the  solar  and  infrared  fields  has  been  made.  The  goal  of  this  study 
has  been  to  provide  recommendations  to  the  SWOE  modeling  community  on  the  appropriate  techniques 
to  use  for  calculating  the  solar  and  infrared  inputs  to  the  background  energy  budgets. 
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BRIEFING  OUTLINE 


•  Background  and  Purpose  of  the  Research 

•  Solar  and  Infrared  Modeling  Requirements  of  the  BTI/SWOE  Program 

•  Review  of  Representative  Techniques 

•  Results 
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Summary  and  Recommendations 
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BACKGROUND  AND  PURPOSE 
OF  THE  RESEARCH 


•  BACKGROUND 

-  A  Goal  of  the  Balanced  Technology  Initiative  (BTI)  on  Smart  Weapons 
Operability  Ehancement  (SWOE)  is  to  Model  the  Radiant  Field  From 
Complex  Natural  Backgrounds. 

-  One  Must  be  Able  to  Model  the  Thermal  Structure  of  the  Background. 

-  Two  of  the  Inputs  to  the  Energy  Budget  are  Solar  and  Infrared 
Radiation. 


•  PURPOSE  OF  THE  RESEARCH 

-  Review  Approaches  to  Calculate  the  Solar  and  Infrared  Fields  and 
Provide  Recommendations  to  the  BTI/SWOE  Modeling  Community  on 
Appropriate  Techniques. 


RADIATION  SCENE  STRUCTURE 


SOLAR  AND  INFRARED  MODELING 
REQUIREMENTS  FOR  BTI/SWOE 


•  FLUX  INPUTS  ARE  REQUIRED  TO: 

-  Thermally  Load  Scenes  Prior  to  Time  of  Scene  Simulation 

-  Calculate  Radiant  Components  Along  Arbitrary  Lines-of-Sight  At  Time 


of  Scene  Simulation 

•  MODELS  MUST: 

-  Valid  for  Wide  Range  of  Environmental  Conditions 

-  Account  for  Cloud  Reflections,  Emissions,  and  Shadowing 

-  Allow  Variable  Atmospheric  Composition  (H2O,  Aerosols,...) 

-  Include  Multiple  Reflections  and  Emissions  from  Scene  Elements 

•  INPUTS  MUST: 

-  Simple  Enough  for  the  SWOE  User  to  Understand 

-  Easily  Obtainable  (i.e..  Surface  Weather  Observations) 


SOLAR  MODELS 
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COMPARISON  OF  SOLAR  TREATMENTS 


SPAF.rA  INC. 


•  Techniques  to  be  Compared: 

-  Solar  Treatment  in  the  Terrain  Surface  Temperature  Model  (TSTM) 

-  ILUMA  (EOSAEL  87) 

-  LOV\7TRAN  7 


•  Results  will  be  Compared  for  Clear  and  Cloudy  Sky  Conditions 


SOLAR  TREATMENT  IN  TSTM 
GENERAL  BACKGROUND 


•  Based  on  Early  Techniques  Used  by  General  Circulation  Models 

•  Assumes  Rayleigh  Scattering  Dominates  Below  0.9  ^m  and 
H2O  Absorption  Dominates  Above 

•  Does  Not  Account  for  Aerosol  Effects 

•  Includes  Adustments  for  Cloud  Cover  and  Slope  of  Receiving  Surface 

•  Cloud  Adjustment  Tied  to  Cloud  Type 

•  Does  not  Separate  Direct  and  Diffuse  Components 


SOLAR  TREATMENT  IN  ILUMA 
GENERAL  BACKGROUND 


SPARTA  INC. 

•  Based  on  Shapiro  Model  (1982)  for  Calculating  Solar  Fluxes  From 
Standard  Surface  Meteorological  Observations 

•  Input  Parameters  Tied  to; 

-  Location 

-  Time  of  Day 

-  Type  of  Surface 

-  State  of  Surface  (Wet,  Dry,  Frozen,...) 

-  State  of  Weather  (Precip,  Clouds,  Presence  of  Obscurants,...) 

•  Not  Tied  to  Temperature  or  RH 

•  No  Accounting  for  Aerosols 

•  Does  Not  Separate  Direct  and  Diffuse  Components 

•  Model  Evaluated  for  Limited  Locations  and  Seasons 


SPARTA  INC. 


COMPARISON  OF 
SOLAR  MODELING  FEATURES 


FEATURE 

TSTM 

ILUMA 

LOWTRAN  7 

Wavelength  Dependent 

No 

No 

Yes 

Separate  Direct  and 
Diffuse  Terms 

No 

No 

Yes 

Adjustments  for 

Clouds 

Yes 

(Empirical) 

Yes 

(Implicit) 

Yes 

(Explict) 

Cloud  Layering 

No 

Yes 

No 

Calculates  Directional 
Radiances 

No 

No 

Yes 

Includes  Weather  Effects 

No 

Yes 

(Implicitly) 

Yes 

(Explicitly) 

RESULTS  OF  COMPARISONS 
SOLAR  FLUXES 


CLEAR  SKIES.  10  SEPT  1982 


TIME  (LST) 
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RESULTS  OF  COMPARISONS 
SOLAR  FLUXES 


CLOUDY  SKIES,  7  SEPT  1982 
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SURFACE  TEMPERATURE  (C) t 


LVOtk  two 


RESULTS  OF  COMPARISONS 
IMPACT  OF  SOLAR  FLUX  DIFFERENCES 
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CLEAR  SKIES,  10  SEPT  1982 


TIME  (LST) 


INFRARED  MODELS 


70 


060*00 


COMPARISON  OF  INFRARED  TREATMENTS 


FEATURE 

TSTM 

Martin  & 
Berdahi 

Schmetz 
et  at. 

LOWTRAN  7 

Clear  Sky  Emissivity 

Tied  To 

Air  Temp.  RH 

Dew  Point  Temp 

Air  Temp 
(Idso  &  Jackson) 

Model 

Atmosphere 

Adjustments  for 

Clouds 

Yes 

(Empirical) 

Yes 

(Empirical) 

Yes 

(Empirical) 

Yes 

(Explicit) 

Cloud  Layering 

No 

No 

No 

No 

Calculates  Directional 
Radiances 

No 

No 

No 

Yes 

Wavelength  Dependent 

No 

No 

No 

Yes 

Includes  Weather  Effects 

No 

No 

No 

Yes 

(Explicitly) 

RESULTS  OF  COMPARISONS 
IR  FLUXES 


CLEAR  SKIES,  10  SEPT  1982 


TIME  (LST) 


RESULTS  OF  COMPARISONS 
IR  FLUXES  -  CLOUDY  SKIES 


Cumulus  Clouds 

(W/nf )  Summer  Midlatitude  Winter 


SUMMARY  AND 
CONCLUSIONS 


^  SUMMARY  OF  MODELING  REQUIREMENTS 
\  FOR  BTI/SWOE 

SPARTA  INC. 


AREA 

REQUIREMENT 

Energy  Budget  Calculations 

Direct  &  Diffuse  Solar  Terms 

Account  for  Reflections  Off  Scene  Elements 

Account  for  Absorption  &  Scattering  by 
Atmospheric  Gases,  Aerosols,  &  Clouds 

Cloud  Contributions  Including  Shadowing 

Radiant  Field  Calculations 

Wavelength  Dependence 

Directional  Scattering,  Reflection,  and 

Emission  Along  Arbitrarty  Lines-of  Sight 

Environmental 

Valid  for  All  SWOE  Scenarios 

Data 

Obtainable  From  Routine  Surface  Weather 
Observations 

SUMMARY  OF  SOLAR  MODELS  STUDIED 
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MODELING 

REQUIREMENT 

TSTM 

ILUMA 

LOWTRAN  7 

ENERGY  BUDGETS 

Wavelength  Dependent 

No 

No 

Yes 

Separate  Direct  and 

Diffuse  Terms 

No 

No 

Yes 

Adjustments  for  Clouds 

Yes 

Yes 

Yes 

Cloud  Layering 

No 

Yes 

No 

Cloud  Shadowing 

No 

No 

No 

includes  Aerosol  Effects 

No 

No 

Yes 

RADIANT  FIELDS 

Calculates  Directional 
Radiances 

No 

No 

Yes 

Calculates  for  Specific 
Wavebands 

No 

No 

Yes 

ENVIRONMENTAL 

Valid  for  SWOE  Locations 

?? 

?? 

?? 

DATA 

Uses  Conventional 

Surface  Weather  Data 

Yes 

Yes 

Yes 

Uses  Radiosonde  Data 

No 

No 

Yes 

SUMMARY  OF  INFRARED  MODELS  STUDIED 


MODEUNG 

Martin  & 

Schmetz 

REQUIREMENT 

TSTM 

Berdahl 

et  al. 

LOWTRAN  7 

ENERGY  BUDGETS 

Wavelength  Dependent 

No 

No 

No 

Yes 

Adjustments  for 

Clouds 

Yes 

Yes 

Yes 

Yes 

Cloud  Layering 

No 

No 

No 

No 

Includes  Weather  Effects 

No 

No 

No 

Yes 

RADIANT  FIELDS 

Calculates  Directional 
Radiances 

No 

No 

No 

Yes 

Calculates  for  Specific 
Wavebands 

No 

No 

No 

Yes 

ENVIRONMENTAL 

Valid  for  SWOE 

Locations 

?? 

?? 

?? 

?? 

DATA 

Uses  Conventional 
Surface  Weather  Data 

Yes 

Yes 

Yes 

Yes 

Uses  Radiosonde  Data 

No 

No 

No 

Yes 

CONCLUSIONS  AND 
RECOMMENDATIONS 


CONCLUSIONS 

•  Solar  and  Infrared  Models  Sampled  From  the  Research  Community  do 
not  Satisfy  BTI/SWOE  Requirements 


RECOMMENDATIONS 

•  An  Approach  Similar  to  LOWTRAN  7  Should  be  Used  to  Provide  Solar 
and  IR  Inputs  to  BTI/SWOE 

•  Modifications  and  Enhancements  Required,  but  Within  the  Framework  of 
LOWTRAN  Concept 


Validation  Studies  With  Field  Data  Mandatory 


CONCLUSIONS  AND 
RECOMMENDATIONS 


CONCLUSIONS 

•  Solar  and  Infrared  Models  Sampled  From  the  Research  Community  do 
not  Satisfy  BTI/SWOE  Requirements 


RECOMMENDATIONS 

•  An  Approach  Similar  to  LOWTRAN  7  Should  be  Used  to  Provide  Solar 
and  IR  Inputs  to  BTI/SWOE 

•  Modifications  and  Enhancements  Required,  but  Within  the  Framework  of 
LOWTRAN  Concept 


Validation  Studies  With  Field  Data  Mandatory 


PCTRAN  7  —  AN  IMPLEMENTATION  OF  THE  GL’s 
LOWTRAN  7  MODEL  FOR  THE  PERSONAL  COMPUTER 

J.  Schioedcr 

Omar  Corporation,  129  Univcreity  Road,  Brookline,  MA  02146 

PCTRAN  7  is  an  implementation  of  the  Geophysics  Laboratory’s  LOWTRAN  7  model,  and  as¬ 
sociated  software,  for  the  IBM,  and  compatible,  family  of  personal  computers.  The  package  contains 
software  for  setting  inputs,  help  capability  for  LOWTRAN  parameters,  viewing  of  output  files,  screen 
graphics  and  hard  copy  graphics.  The  software  has  been  validated  by  the  GL  under  a  cooperative  IR  &D 
agreement  with  Ontar. 

This  paper  will  describe  PCTRAN  7,  Version  2,  and  demonstrate  the  capabilities  of  the  package. 


PCTRAN  7  [c]:  An  Implementation  of  the  GL’s 
LOWTRAN  7  Model  for  the  Personal  Computer 


ANNUAL  REVIEW  CONFERENCE  ON  ATMOSPHERIC  TRANSMISSION  MODELS 
Geophysics  Laboratory,  Hanscom  AFB,  MA 
5  June  1990 

Paul  V.  Noah,  and  John  Schroeder 

Ontar  Corporation 

129  University  Road 

Brookiine,  MA  021 46  *4532 

Tel:  617-739-6607  FAX:  617-277-2374 
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Cooperative  R  &  D  Agreement 


With  Geophysics  Laboratory  -  Hanscom  AFB,  MA 
September  1988 

PC  Implementation  of  a  Software  Package-  LOWTRAN  7  -  PCTRAN  7  [c] 


%WUU  IlflABB 


PCTRAN  7  [c]  *  Version  2 


PC  Version  of  the  AFGL  LOWTRAN  7  Atmospheric  Radiance  &  Transmission  Code 
Complete  Implementation  of  the  LOWTRAN  7  Code  -  7.37 
Interactive  User  Input  Software 
Help  Screens  for  All  Input  Variables 
Screen  and  Hard  Copy  Graphics  Output 
Tabular  Output  in  ASCII  Format 
Batch  Processing  Input  Software 
Plotting  from  Different  LOWTRAN  Calculations 

^CERTIFIED  for  ACCURACY  by  the  GEOPHYSICS  LABORATORY 


Hardware  and  Software  Requirements 


Personal  Computer  -  XT,  AT,  80386,  80486  (Compatible,  Clone) 

1.2  Mbyte  Diskette  Drive,  Hard  Disk 
640  Kbytes  of  Memory 

CGA,  EGA,  or  VGA  Graphics  Board  and  Monitor  -  for  Screen  Plots 

Printer  -  for  Hard  Copy 

Numeric  Co-processor  Highly  Recommended 

i 
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LITERATURE  AVAILABLE 


SOFTWARE  DEMONSTRATION 


Ontar  Corporation 
129  University  Road 
Brookiine,  MA  02146 -4532 

Tei:  617-739-6607  FAX:  617-277-2374  Buiietin  Board:  617-277-6299 


ONTAR’s  LOWTRAN  Program  Suite  -  Version  7.2 


a.  LOWTRAN  input  -  LOWIN 

b.  Execute  LOWTRAN 

c-  LOWTRAN  Piotting  -  LOWPLT 

d.  Printer  Plotting  -  PRTPLT 

e.  View  LOWTRAN  Output  -  VIEWOUT 

f.  View  FILE7  -  VIEWOUT 

g.  View  FILES  -  VIEWOUT 

h.  LOWFIL  input  -  LFILIN 

i-  Execute  LOWFIL 

).  View  LOWFIL  Output -VIEWOUT 

k.  Execute  LOWSCAN 

l.  LOWSCAN  Plotting  -  LOWPLT 

m.  View  FILE9~  VIEWOUT 

n.  Multiple  LOWTRAN  Plots  -  LOWMPIN  &  LOWPLT 


o.  Return  to  DOS 


Input  Function: 


PC-TRAN7  Batch  Mode  Manager. 


ESC  -  Quit  LOWIN  (write  LOWIN  and  LOWPLT.DAT) 

F2  -  Edit  current  run. 

F3  -  Edit  next  run. 

F4  -  Edit  previous  run. 

F5  -  Add  new  run  (to  end)  and  go  to  that  run 
F6  -  Delete  current  run. 

F7  -  Go  to  run. 

Database  name  MCASE3 

Number  of  runs  In  this  database  4  Current  run  3 


LOWTRAN7  Card  5 


Initial  Altitude  (km)  1 .500 

Final  Altitude/T angent  Height  (km)  1 0.000 

Initial  Zenith  Angle  (degrees)  .000 

Path  Length  (km)  .000 

Earth  Center  Angle  (degrees)  .000 

Radius  of  Earth  (km)  [.000  -  default]  .000 

Type  of  Path  Short 


Initial  Frequency  2000.000  cm-1  Wavelength  4.000  m 
Final  Frequency  2500.000  cm-1  Wavelength  5.000  m 
Frequency  Increment  (wavenumber)  5.000 


Run  #  2  of  4 


LOWTRAN7 


Cards  3  &  4 


Plot  Type 

Transmittance  in  m 

Type  of  X  Axis 

Linear 

Type  of  Y  Axis 

Linear 

Number  of  Decimal  Digits  for  Y  Axis 

2 

Length  of  X  Axis  (In  inches) 

7.0000 

Beginning  Wavenumber/Wavelength 

4.0000  m 

Ending  Wavenumber/Wavelength 

5.0000  m 

X  Axis  Annotation  Interval 

.2000  m 

X  -  Number  of  Minor  Ticks  /  Division 

5 

Length  of  Y  Axis  (in  inches) 

6.0000 

Autoscale  Y  Axis 

No 

Minimum  Transmittance/Radiance 

.OOE+00 

Maximum  Transmittance/Radiance 

1.00E+00 

Y  Axis  Annotation  Interval 

2.00E-01 

Y  -  Number  of  Minor  Ticks  /  Division 

5 

Plot  Grids  (Graph  Paper) 

Coarse  Grid 

Run  #  2  of  4  LOWPLT  Scaling  Plot  #  1 
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Multiple  Run  Plotting  Inputs 


Filename 

Run 

Plot  Mode 

Title 

CASE-A 

1 

Transmittance 

1976  U  S  STANDARD 

CASE-A 

2 

Transmittance 

SUBARCTIC  WINTER 

CASE-A 

3 

Radiance  /w  Scattering 

SUBARCTIC  SUMMER 

CASE-A 

4 

Radiance 

MIDLATITUDE  SUMMER 

CASE-B 

1 

Transmittance 

1976  U  S  STANDARD 

CASE-B 

2 

Radiance 

TROPICAL  MODEL 

CASE-B 

3 

Radiance  /w  Scattering 

MIDLATITUDE  SUMMER 

CASE-B 

4 

Transmittance 

SUBARCTIC  SUMMER 

CASE-C 

1 

Transmittance 

1976  US  STANDARD 

CASE-C 

2 

Radiance 

TROPICAL  MODEL 

CASE-C 

3 

Radiance  /w  Scattering 

New  Model  Atmosphere 

CASE-C 

4 

Transmittance 

Met  Data  (Hor  Path) 

LOWMPIN  Select  Runs 


ROD  I RNCE  110'^  WRTTS/CM2- 


UNVtLLNGTH  (M  I  CRCLMG  1  hfi) 

.  1  y .  13  9.7  9.  G  9.5  9.4  9.9  9.2 


WAVENUMBER  lCM-1) 


UJ 

CJ 

2: 

CE 


cn 


cr 

cn 


MULTIPLE  PLOT  TEST  CFISE 
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SPECTRAL  MODELING  OF  OFF-AXIS  LEAKAGE 
RADIANCE  USING  THE  MODTRAN  CODE 
N.  Grossbard 

Boston  College,  885  Centre  Street,  Newton,  MA  02159 
D.R.  Smith 

Geophysics  Laboratory,  Hanscom  AFB,  MA  01731 

MODTRAN  has  been  used  together  with  the  Degges  High  Altitude  Infrared  Radiance  Model  for 
the  purpose  of  modeling  the  off-axis  leakage  spectra  obtained  from  the  HIRIS  and  SPIRE  nocketbome 
atmospheric  experiments.  Off-axis  leakage  spectra  for  numerous  angular  shells  around  tltc  sensor 
FOV  were  calculated  using  nominal  or  pre-flight  estimates  of  the  telescope’s  off-axis  rejection  (OAR) 
pcrformarKC.  Each  leakage  contribution  was  weighted  to  obtain  a  best  least  squared  fit  to  the  actual 
flight  data.  The  resulting  weighting  functions  were  then  applied  as  correction  factors  to  the  initial 
estimates  of  the  off-axis  rejection  performance  to  obtain  estimates  of  the  inflight  OAR  performance  of 
the  telescope.  The  computed  spectra  arc  in  good  agreement  with  the  flight  data  in  each  case  and  the 
resulting  OAR  performance  curves,  though  considerably  worse  than  initial  estimates  based  on  laboratory 
measurements,  are  reasonable  for  inflight  conditions  and  in  good  agreement  with  other  inflight  estimates 
of  OAR  performance. 


Spectral  Modeliiig  of  Off-Axis  Leakage  Radiance  Using  the  MODTRAN  Code 
N.  Grossband  (Boston  College)  and  D.R.  Smith  (Geophysics  Laboratory) 
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CLOUD  OPACITY  RETRIEVAL  USING  LOWTRAN 
AND  THE  STAMNES  SCATTERING  MODEL 
B.L.  Lindner  and  R.G.  Isaacs 

Almospheric  and  Environmental  Research,  fnc.,  840  Memorial  Drive,  Cambridge,  MA,  02139 

The  AFGL  LOWTRAN  atmospheric  transmittance  model  and  an  efficient  Mie  theory  algorithm 
have  been  combined  with  the  Stamnes  discrete  ordinate  method  multiple-scattering  model,  and  used  to 
simulate  multispcctral  sensor  intensities  for  a  variety  of  cloud  properties,  surface  types,  sun  and  sensor 
geometries,  and  background  atmospheres,  from  visible  to  infrared  wavelengths.  These  intensities  are 
stored  in  a  lookup  table,  and  a  minimization  procedure  has  been  developed  to  most  efficiently  and 
accurately  match  observed  intensities  with  those  in  the  lookup  table,  and  hence  retrieve  cloud  properties. 
The  minimization  procedure  has  been  used  to  retrieve  cloud  optical  depth  from  LANDSAT  TM  data. 


COMBINE  MODELS; 


u 

Q 

O 

U 

O 

z 

M 

cd 

u 


u  < 

U  O 
O  CO 

X 

CO  u 


td 

a 

o 

o 


o 


CO 


0^ 

UJ 

N 

<0 

W 

X 

y-N 

0 

CO 

Id 

hH 

(d 

0 

H 

cd 

2 

UJ 

< 

z 

CXi 

w 

Oi 

,4 

CO 

2 

“ 

0 

0 

CO 

z 

£ 

CO 

< 

H 

0 

LO 

< 

Id 

z 

3 

H 

Cd 

Cd 

0 

CO 

*— * 

> 

< 

p 

UJ 


CO 

UJ 

(d 

X 

UJ 

l-H 

0 

Cd 

H 

UJ 

0^ 

Id 

0 

p 

ui 

CLi 

-  (Cl 

oi 

Cd 

H  0 

H 

0 

a 

0 

CO 

w 

2 

E  fd 

Cd 

z 

z 

Cl 

U 

UJ 

Pd 

A  Q 

< 

CO 

CO 

u 

J  P 

u< 

■ 

0 

►  0 

z 

X 

5 

5 

H 

c  CJ 

CO 

CO 

< 

< 

H 

< 

Q 

H 

ct: 

(d 


o 

< 

a, 

o 

Q 

O 

hJ 

o 


u 

a 

o 

2: 

o 

z 

KH 

o: 


CO 

o 


j= 

o 

V4 

<1) 

w 

Cci 


Ld 

0 

H 

u 

H 

Cli 

c 

< 

(U 

U 

*0 

CO 

c 

0 

CO 

V4 

Cd 

cd 

z 

Cd 

> 

£ 

c 

<. 

Q 

Cd 

H 

2 

CO 

HH 

T? 

X 

c 

Cd 

(3 

Z 

X 

H 

CJ 

cQ 

Q 

V4 

0^ 

X 

a. 

(0 

o 

e 


•if, 


radiance  calcuuted  bt  solving. 


(p»nuv3uo3)  oHissaDoaa  viva 


48 


AtmospherLc  Vertical.  Opacity  at  Solccti-.i 
Wavelengths  for  Midlatitude  Summer  Condi tior. 


e«r«  Soil.  SoUr  ZonUh  Anglo  of  61* 


CALCULATED  CLOUD  EDGE  SWIR  RADIANCE  GRADIENTS  AS  VIEWED  BY  SATELLITE 

L.L.  Sm’th 

Grumman  Corporate  Research  Center,  Bethpage,  NY  11714 

As  surveillance  systems  become  more  capable  and  more  sophisticated,  knowledge  of  the  effects  of 
clouds  as  background  clutter  becomes  an  irKteasing  stressing  factor  on  performance.  Cloud  edge  radiance 
gradients  in  the  SWIR  arc  calculated  using  the  LOWTRAN  7  code  for  multiple  satellite -cloud-sun  viewing 
geometries.  As  expected,  the  gradient,  defined  as  the  change  in  cloud  radiance  per  degree  in  the  detector 
focal  plane,  increases  as  the  viewing  position  approaches  the  earth’s  limb.  This  is  principally  caused  by 
the  foreshortening  of  the  cloud  boundary  as  the  limb  is  approached  with  a  small  effect  due  to  an  increase 
in  cloud  brightness.  The  gradient  peaks  in  both  the  forward  scattering  and  backscattering  directions  arc 
due  to  the  phase  angle  effect  of  the  scattering  particles. 


Calculated  Cloud  Edge  SWIR  Radiance 
Gradients  as  Viewed  by  Satellite 


Lewis  L.  Smith 

Grumman  Corporate  Research  Center 
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Scattering  Angle  Definition 


SCATTERING  ANGLE  =  COS  'fsiN  SIN  (Zs^^)  COS  -  COS  COS  (Zs^f^] 


Cloud  Edge  Gradient  Model 


•  SWIR  radiances  calculated  using  standard 
cirrus  cloud  model  from  L0WTRAN7  w/multiple 
scattering  &  daytime  conditions 

•  Cloud  thickness  varies  in  step  size 

-  1 1 1  km  in  longitude  (east) 

-  1 1 1  km  in  latitude  (south) 

•  Maritime  aerosols 

•  Sub-satellite  point  on  equator 


Cloud  Edge  Gradient  Model  (Cont'd) 


.  Model  atmospheres 

-  0-25°  tropical 

-  25-55°  midlatitude  summer 

-  55-70°  subarctic  summer 

.  Cirrus  cloud  base  altitude 

-  0-25°  1 1  km 

-  25-55°  1 0  km 

-  55-70°  7  km 

•  Cloud  edge  gradients  calculated  from 
changing  radiance  in  satellite  detector  focal  plane 

.  All  gradients  normalized  to  cloud  at  sub-satellite  point 


Cloud  Gradients  -  South 
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Cloud  Gradients  -  East 
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Cloud  Gradients  -  East 
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Conclusions 


•  Cloud  edge  SWIR  gradient  increases  as  viewing 
position  approaches  earth's  limb  -  caused  by 
foreshortening 

•  Gradient  increases  as  cloud  altitude  increases  due 
to  decreased  atmospheric  absorption 

•  Gradient  peaks  in  both  forward  and  backscattered 
directions  due  to  phase  angle  effect  of  scattering 
particles 


THE  DEPARTMENT  OF  ENERGY  INITIATIVE  ON  ATMOSPHERIC  RADIATION 
MEASUREMENTS  (ARM):  A  STUDY  OF  RADIATIVE  FORCING  AND  FEEDBACKS 

R.G.  EUingson 

Department  of  Meteorology,  University  of  Maryland,  College  Park,  MD  20742 

G.M.  Stokes 

Applied  Physics  Center,  Pacific  Northwest  Laboratories,  Richlarxi,  WA  99352 

A.Patrinos 

Division  of  Atmospheric  and  Qimate  Research,  US  Department  of  Energy,  Washington,  DC  20545 

As  a  key  component  of  the  strategy  to  address  global  climate,  the  Department  of  Eneigy  has  launched 
the  Atmospheric  Radiation  Measurements  (ARM)  initiative.  The  objectives  of  the  ARM  Program  are  to 
provide  detailed  measurements  of  radiative  effects  in  the  atmosphere,  and  to  provide  parametcrizations 
of  these  effects  for  use  in  atmospheric  models.  Particular  emphasis  will  be  on  Uiose  effects  associated 
with  clouds  and  grccnhou.se  gases.  The  effort  will  support  the  continued  and  rapid  improvement  of 
GCM  predictive  capability.  The  presentation  will  summarize  the  science  context,  program  requirements, 
measurement  strategy,  experimental  approach,  scientific  management  and  site  selection. 


State  of  the  Art 
in  Radiation  Modeis 

ICRCCM  has  shown  a  wide  range  of 
disagreement  among  radiation  models  used  in 
climate  studies 

There  are  no  calibrated  techniques  for  predicting 
and  calculating  the  radiative  effects  of  clouds  for 
either  homogeneous  or  broken  conditions 


Radiation  and  Climate  - 
Important  Facts 

Radiation  is  a  quantity  to  which  Earth's 
climate  is  VERY  sensitive  - 

1  %  changes  are  important 


Atmospheric  Radiation 
Measurement  Program 


A  Study  of  Radiative  Forcing 
and  Feedbacks 


General  Goals  of  ARM: 


Improve  the  Performance  of 
General  CirculationModels 
of  the  Atmosphere 
as  Tools  for  Predicting 
Global  and  Regional  Change 


Specific  Goals  of  ARM 


Improve  the  Tk*eatment  of 
Radiative  transfer  in  GCMs  Under 
Clear  Sky,  General  Overcast,  and 
Broken  Cloud  Conditions 

Improve  the  Parameterization 
of  the  Properties  and  Formation 
of  Clouds  in  GCMs 


ARM  Design  Strategy 

Succcssi  \  ’c  A  pproximation 

Scientific  Issues 
Experiment  Definition 
Instrument,  Model,  Site 
Selection  (Process) 

Data  System  Requirements 


Experiment  Definition 

Key  Organizing  Concept 

Model  to  be  Tested  (hypothesis) 
Input  Data 
Test  Data 


Experimental  Context 


Hypothesis-Oriented  Radiation 
Science  Experiment 

Clouds  and  Radiation  Testbed 


Model  Definition 


Numerical  (computational) 
representation  of  the 
understanding  of 
a  physical  process 


Design  Process  Goals 


Rapid  Deployment 
Long-term  Flexibility 
Strong  Coupling  to  Science 


CART  Data  System: 
Design  Activities 


Data  System  Design: 
Objective 


A  data  system  that  can  provide 
streams  of  data  (measurements  and 
predictions)  of  known  quality  that 
can  be  compared  to  test  the  predictive 
power  of  radiometric  and  other 
models  selected  to  be  a  part  of  the 
ARM  program 

i 

i _ 


Requirements  Definition 


•  Formal  top-down  design  sessions  to  include 
representatives  of  Science  Team  and 
instrument  developers 

•  Design  reviews  and  workshops  with 
Science  Team  and  interagency  working 
groups 

•  Formal  design  review  with  outside  panel 


Summary  of  Design  Process 

•  Design  driven  by  scientific  requirements 

•  Focus  on  providing  streams  of 
measurements  and  predictions  for 
intercomparison 

•  Creation  of  a  process  that  supports 
evolution  of  ARM  scientific  activities 


THE  HITRAN  MOLECULAR  DATABASE  IN  1990 
Lt.  S.  Shannon  and  L.S.  Rothnian 
Geophysics  Laboratory/OPI,  Hanscom  AFB.  MA  01731-5000 

The  spectroscopic  molecular  database,  HITRAN.  is  the  DoD  and  international  standard  compilation  of 
absorption  parameters  that  enable  the  calculation  of  atmospheric  spectral  simulations  from  the  microwave 
through  the  visible.  HITRAN  has  been  periodically  improved  and  released;  the  current  edition  being  1986 
(Applied  Optics,  26,  4058  (1987)).  This  present  edition  contaias  over  a  third  of  a  million  transitions  for 
some  28  species  and  their  significant  atmospheric  isotopic  variants.  The  new  edition  will  have  additional 
species  and  will  contain  many  more  transitions.  Furthermore,  there  will  be  appended  new  cross-sections 
for  heavy  molecular  species,  with  bands  at  several  reprc.scntativc  temperatures.  This  talk  will  summarize 
some  of  the  major  updates  and  modifications  that  will  be  available  on  HITRAN  1990. 


The  HITRAN  Molecular  Database  in  1990 


Lt  Scott  Shoi  inon  &  Laurence  S.  Rothman 
GL/OPl 

Hanscom  AFB.  MA  01731-5000 


OUTLINE 

0  Major  Molecular  Parameter  Updates 
®  SELECT:  Database  Interface 
(D  Conclusions 


Major  Molecular  Parameter  Updates 

O  HjO 
O  CO2 
O  O3 
O  CH4 

O  CO,  HNO3,  HCI,  HF.  HI.  HBr... 


O  Cross-sections  (CIONO2,  CFC’s,  N2O5...) 
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Carbon  Dioxide  (CO2) 

O  New  Line  Positions:  New  fit  of  observed  high  resolution  data 
by  Hawkins  &  Rothman.  New  data  includes  very  high  vibrational  4.3- 
obsen/ations  of  Bailly  et  al.  (Orsay)  and  the  15-tim  observations 
of  Esplin  et  at.  (GL) 

©  New  intensities:  New  intensity  observations  of  Dana  et  cd. 
(Paris)  and  Johns  (NRC).  Calculations  by  Wattson  &  Rothman  for 
unobserved  bands,  include  Herman-Wallis  coefficients. 


®  Revised  Halfwidths:  Re-evaluation  of  all  halfwidth  observa¬ 
tions. 


Ozone  (O3) 

O  New  data  for  V2  and  2v2  -  V2  from  Pickett  (JPL)  et  aL 

©  Umpteenth  revamping  of  10-|im  region.  This  time  includes 

many  new  bands  and  hot  bands. 

©  Isotopic  bands  of  10-|im  region  from  FCP  &  Rinsland  (NASA 
Langley) 

O  Many  new  combination  bands  from  Goldman  (Denver) 
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Methane  (CH4) 

O  Major  improvement  by  Brown  (JPL) 

©  Updated  mono-deuterated  methane  bands  by  Brown. 
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Other  Molecular  Species 

O  CO  (carbon  monoxide):  Update,  primarily  for  HITEMP  (Tipping) 

O  NO  (nitric  oxide):  Update  of  fundamental  (Ballard) 

O  HNO3  (nitric  acid):  addition  of  new  bands 

(fundamentals  and  combinations,  410  to  1400  cm-^) 
(Goldman,  Maki,  Rinsland,  Perrin,  etal.) 

O  HCI,  HF,  HI,  HBr  (hydrogen  halides): 

Update,  primarily  for  HITEMP  (Tipping,  etal.) 

O  H2O2  (hydrogen  peroxide);  Vg  band,  1 38-1500  cm  ’  (Hillman) 

O  C2H2  (acetylene):  vg  band,  1192-1470  cm-’  (Rinsland) 

O  C2H6  (ethane),  vg  and  V7  bands,  0-833  and  2973-3000  cm  ’ 
(Blass:  Goldman,  Dang-Nhu,  etal.) 

O  COF2,  SFe...:  new  species  on  compilation  (Goldman). 
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Cross-sections  (Supplemental  files) 

CFC’s  now  at  six  temperatures:  203,  213,  233,  253,  273,  and  293K. 
N2O5  at  233,  253,  273,  and  293K.  CIONO2  at  21 3  and  296K. 
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Cross-sections  (continued) 
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SELECT  OPTIONS 

Input: 

(X)  V,  and  V2  (Initial  and  final  frequencies  of  selection) 

®  Molecule 
<S)  Isotope 

®  v',  v"  (upper  and  lower  "global"  quanta) 

®  Sent  (Intensity  cutoff) 

Output: 

O  Batch  file  for  input 
O  Files  for  input  to  subsequent  programs 
(direct  image;  '82  format;  user  defined) 

Q  Hard  copy  listing 

(codes  converted  to  spectroscopic  and  chemical  notations; 
80  or  132  columns  printers) 


SUMMARY 


•k  Improvements  notable  for  remote  sensing,  laser  transmission, 
climate  modeling,  use  of  heavy  species... 

★  New  Edition  of  HITRAN  to  appear  in  fall  1990. 

k  Users  on  mailing  list  will  be  notified  through  newsletter. 

k  HITEMP  due  snortly  after  HITRAN. 

SELECT  faster  and  more  flexible. 


HITRAN  Database  Distribution 
Tape;  9  irack.  ASCII.  1600  BP! 

Ciimafological  Services  Section 
National  Climatic  Data  Center  of  NOAA 
Federal  Building 
Asheville,  NC  28801 
phone:  (704)  259-0682 

or 

Or.  Laurence  S.  Rothman 
AFGL/OPI 

Honscom  AFB,  WA  01731-5000 
phon^:  (617)  577-2856 


LINE  COUPLING  CALCULATIONS  FOR  INFRARED 
BANDS  OF  CARBON  DIOXIDE  FOR  FASCOD3 
M.L.  Hoke,  F.X.  Kneizys,  S.A.  Qough*,  J.H.  Chclwynd 
Geophysics  Laboratory/OP,  Hanscom  AFB,  MA  01731 

Second  order  perturbation  theory  line  coupling  parameters  (y’s  and  g’s)  have  been  calculated  for  the 
six  most  important  Q-branches  of  carbon  dioxide  in  the  fifteen  micron  spectral  region;  those  at  618,  667 
(fundamental),  720,  741  and  791  cm  *  of  the  most  abundant  isotope  016-C12-016  and  at  721  cm  *  of 
the  less  abundant  isotope  016-C12-O17.  Line  coupling  parameters  have  also  been  calculated  for  the 
P  and  R  branches  of  the  4.3  micron  fundamental  band.  For  each  band  the  temperature  dependence  of 
the  coefficients  is  accounted  for  by  interpolation  of  a  set  of  four  pairs  of  coupling  parameters  (y,g)  at 
the  temperatures  200,  250,  296  and  340  K.  These  line  coupling  parameters  are  part  of  the  new  version 
(version  3)  of  the  GL  high  resolution  line-by-line  computer  code  FASCODE.  Details  of  the  calculations, 
including  the  appropriateness  of  perturbation  theory  and  the  accuracy  of  the  temperature  interpolation 
scheme,  will  be  discussed. 

*Now  at  Atmospheric  and  Environmental  Research,  Inc.,  840  Memorial  Drive,  Cambridge,  MA 
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FASCOD3:  AN  UPDATE  (WITH  NLTE  EMPHASIS) 

G.P.  Anderson,  F.X.  Kneizys,  E.P.  ShetUc*,  J.H.  Chetwynd,  L.W.  Abrcu,  M.L.  Hoke 
Geophysics  Laboratory/OP,  Hanscom  AFB,  MA,  01731 
S.A.  Clough,  R.D.  Worshani 

Atmospherie  and  Environmental  Research,  Inc.,  840  Memorial  Drive,  Cambridge,  MA,  02139 

FASCOD3,  a  line-by-line  atmospheric  radiance-transmittance  code,  is  currently  scheduled  for  a  late 
summer  release.  As  with  FASCOD2,  the  program  is  applicable  to  spectral  regions  from  the  microwave  to 
the  middle  ultraviolet,  employing  standard  spectroscopic  parameters  supplied  from  external  line  aliases. 
New  FASCOD3  capabilities  include  multiple  scattering  of  thermal  radiation,  CO2  and  O2  temperature- 
dependent  line  coupling,  UV  diffuse  absorption  (O2  and  O3),  improved  weighting  functions,  and  enhanced 
■lon-local  thermodynamic  equilibrium  (NLTE)  calculations.  NLTE  input  requirements  include  temperature 
or  population  profiles  for  the  specified  excited  vibrational  slates.  Any  auxiliary  NLTE  line  positions, 
halfwidths,  strengths  and  vibrational-rotational  assignments  must  also  be  provided. 


*Now  at  NRL. 


I'’ASC'()D.<  An  Updiilc 


u 

£1  "ZL 
^  . 

s-  i  < 

c  ^  -e 
..  ■  J 
x:  .3 
^  -  A 

•  "Tj  jC 
c  C  Q. 
$  - 


-S  a: 


$  c 

.  D 

Q  E 
a:  o 


t/5 

"O 


FASCODE:  Fast  Ataospharlc  Slqnatuto  Code 


o 

o 


& 

•ri 

o 

•H 

u 

A 

a 

o 

a 

< 


i4 


n 

A 

ICO 

Q 

iQ 

41  Vi 

O 

a 

c 

U 

<a 

•H  c 

3 

m 

o 

D 

u. 

S 

u 

0  o 

0 

u 

0 

c 

o 

c 

0 

I 

¥> 

» 


Q 

O 

o 


O 

O 

o 


v> 

E 


u 

o 


a  g 


p*  o 
C  Q 
^  O 

o 


s 

z 


55' 


«  >. 
>  «  « 
b-  «  «-- 
•o  O  t7* 

-•  o  o 
oc  X  u 


•a 


-3  <  U 
z  a  ^ 


0 

M 

tt 

U 

V> 

—4 

c 

C 

0 

0 

0 

c 

•H 

■H 

c 

-o 

l4 

u 

0 

>c 

0 

0 

0 

x: 

o 

k4 

x: 

v< 

u 

A> 

VI 

O' 

» 

a 

T3 

o 

0 

c 

m 

a 

C 

0 

•H 

3 

H 

0 

T3 

U 

tt 

A 

« 

C 

c 

c 

0 

O' 

0 

a 

« 

V( 

0 

c 

tt 

U 

tt 

•H 

0 

0 

0 

c 

0 

• 

0 

x: 

a 

AJ 

u 

• 

0 

AJ 

o 

3  — » 

0 

0 

Vt 

0 

A> 

c 

0  y 

9 

0 

A 


«  u  o 

^  c 

x)  x: 

I  U  >4  %4  >  k 

or  •H  0  o  «  0  M 

xi  %4  n  >  « 

t^  ^  C  C  «  O'  c 

C  M  0  O  C  C  •H 

o  n  o  ■ 


a 

TJ 


Q. 

O 


a 

3 

O 

o 


£  tH 

<0  fi  3  a 
C  C  3 
*H  -*4  <0  0 
0  fi  «»  U 

«  C 

4J  «  «  I 
C  > 

O  «  -<  « 
U  ^  V  u 
w  I  u  a> 


o 

V4  *0 

M  C 

«  « 

^  0 

0  0  (M 

*>  O 
>%  lA  O 

•H'4  <a 

C  tt 

H  *0  ^  «Q 

U  C  X3  «  0 

•<H  «  <q  ^  0  ‘*0 

>  Xi  «-4  U  C  <0 

u 

XT 

C  Q. 

«  -4  a 

^  O  4J 

0  «  3 

•p  0  ^  a 

t  IQ  «Q  >-4  r-4  C 

Ul  ^  3  'H 

O  O  IQ 

a  TJ  > 

E  O  O  U 

•*4  E  T3  C 


•O  0 
C  E 


C  t  «  w  o* 


0  C 

a  3 

K  !■« 

£ 

O'  lA  C 

c  o 


■u  c 

C  0 

IQ  TJ 


C  *0  c 

O'  o  u 

O  U  •-«  0 
O 

0  *0  M  Vi 
i4  0  0  0 
—  -D  -H 

•a  VI  «-♦ 

0  •H  <-l 

t3  <Q 
E  O  E 
0  E  VI 


NOlSSIWSNVai 


NOissiriSNvai 


V 

o 

o 


o 

x: 

0 


0 

t  A> 

kl 

A>  U 

3 

-4  0 

Al 

0 

3  x: 

0 

0 

n 

0  u 

x: 

c 

0 

V 

AJ 

O' 

0 

0  y 

c 

•o  o 

o 

« 

VO 

0 

0 

(  r>l 

E 

JC 

VI 

u 

aJ 

Vi 

o 

U  0 

••4 

C 

0 

C  t 

0  *0 

«-» 

D  C 

Ai 

o 

C 

0  0 

c 

o  o 

0 

a 

o 

fi  c 

C  0 

0 

U  *-1 

x: 

3 

O 

0 

3 

to  A-l 

a 

A> 

WO 

-4  0 

«<  u 

V) 

*J  0 

Ai  C 

a 

U.  0  c 

o 

c 

0  U  C 

0  o 

•• 

C  3 

R 

0 

•p4 

3  3  0 

o  -1 

0  O 

VI  0  V( 

AJ 

C 

0  V4 

-4  U  -<-< 

^  A> 

0  X 

O  -4 

< 

o 

C  0 

a  AJ  AJ 

V(  a 

0  3 

AJ  0 

-4 

0  *0 

—  W  -4 

0 

0 

>4  a  c 

a 

0 


V)  t 

0  0 
>  *J 
C  0 
O  E 
O  « 

u 

0  IQ 

c  a 


>  T3  4:  x: 

<Q  a  u 
0  c  -Q 
»  a  X 


V*  0  E 
3  U  0 
Q.  V  XT' 
Vt  *J  to  ^ 
0  3  C 
V)  O  0  « 
IQ  (A  £ 
U  0  IQ 
^  X3  0 
V  O.  W  ^ 
•J)  B  4J  ■'t 
0  IQ  <Q  U. 
H  M  O  ^ 


O  U  H  *J 


O  0 
c 

AJ  "4 

a  VI 

,  ^  3 


•Q  C  C  O  0 

3  IQ  0  0  y  u  ' 


a  —4  V* 

O  0 

to  a 


> 

o 

Xi 

IQ 


fi-i; 


0  -• 

a  a 

to  3 


i4 

O  14 

V»  U  -  . 

«  w  >4  x:  o  3 

0  i3  V}  U  C 

£  T3  » 

0  0  0  aJ 

>4  6  C  C  C  C 

x:  x:  •-«  •-I  o 

a  *J  J  J  iJ  u 

w 

O  U 


-  a 
a 

0  M 

C  O 
O  tn  0 
^  o 
AJ  »-•  0 


U  r-l  C  3  O  l4 

c  0  Vt  ki  <3 


>4tT'0C' 

0  U  C  E  -• 

W  A>  -4 

3  0  u  u  x:  * 

V4  E  0  0  cr 

V4  o  >4  >4  *4  , 

-4  0  «  (Q  0  4 

O  O  *-4  J  IK  ( 


o  (O 
u  H 

0  ij 

<  X  r 


a 

£ 


-10  0  0  0  3 


O 

O' 


u  IQ  t-»  0 

TJ  IQ  V  I' 

•4J  -4  a  0  *-< 

«  -4  r.  V.  c 

0  «  O  0  -••  . 

H  >  o  a  a 


in  fascod:^ 


FASC0D3:  Fast  Atmospheric  Signature  Code 


4J 


c 

o 

V 


t/i 

2 

O 

M 

o 

< 


10 

c 

o 

•H 

u 

c 

2 

O' 

c 

•H 

£ 

cr 

•H 


I 

T3 

V. 


o  o 
u 
a. 
w 


« 


c 

O' 


£ 

a 


£  « 

U 

O  C 

^  « 


o 


a  — 

a 

•  *0 

£  P 

*j  o 


«  e 


§ 

o  a 

^  3 

O' 

c 

X  « 

o 

o 


c 

O  I 
£ 


IQ  «4 

Jii  E 


O* 

C  4; 

«  : 

u  J 


M 

O 

•H 

>  « 

«  *0 

CO  O 
•H  4> 

e  •H 

V  4J 

4H 

•o  « 

c 

«  K 


■* 

e 

•0  -H 

O' 

V) 

M 

4J 

0 

c  V 

c 

a 

« 

X 

3 

0 

«  < 

-H 

0 

b 

M 

a  u 

^4 

A 

0 

« 

0* 

%4 

V) 

V  ^ 

■H 

c 

c 

•w4 

b 

ns 

> 

*w4 

p 

3 

O' 

0: 

>.  *> 

S 

u 

b 

0 

^  0 

« 

O' 

U 

<-<  V 

0 

r* 

c 

U 

4J 

w 

c 

10 

<Q 

0 

.•4 

y 

•w 

X 

C  M 

*0  c 

to 

(> 

T» 

c 

E 

M 

0  V 

Cl  0 

u 

ti 

0 

0 

<0 

irt 

M 

10 

0 

•*-4 

V4 

U 

>N  4J  « 

< 

<3 

i} 

U  V 

Cl  U 

«.* 

u. 

T1 

3 

w 

Cl  u 

u  c 

0 

C 

0 

T1 

u  u 

y  p 

(ft 

N 

Cl 

-4  c 

y  ^ 

c 

0) 

•^4 

3 

0 

a 

0 

*J 

•0 

■D 

u 

tA 

k« 

y 

a 

e 

a 

o 


c  ^  c;  ^ 
o  o  3 
c  u  xi 


•'•*•  •,*"»)/*UCM)  33NVI0VM 


a. 

o 


wavenumber  Ic»-I) 


N'on-loca.  thermodynamic  equilibrium  (NLTE)  -  Radiance 


(J 

o 


fieccnt  Compjrisons/Validation: 


V 

o 

u 

« 

u 

3 

4J 

« 

c 

O' 

10 


« 

u 

AJ 

a> 

o 

£ 


O 


o 

u 

c 

« 


-o 

o 

o; 


V 

C  CO 
<0  GO 

<n 


o  I 
c  m 

lO  CO 


%4  w 

o 


C  £ 
0  x: 


u  u 
«  0 
a  cr 
e 

o  < 

u 


£ 

O' 


^  >  c 
I  «  o 
>1  )  u 
0 

I  U  «  O' 
9  O  JZ  C 

c  u  ^ 
-H  e  O'  T3 
*H  C  3 

^  o  ^ 
'e  c  *H  u 
«  «  c 
x:  >  -H 

«  <0 
^  u  > 

>  -♦  4> 

Tj  E  a 

T3  «  E  <0 

«  c  x: 


Cl.  X3  n  V 
EE  C 

o  a 

(J  M  (-H 


E  T5 

« 

•U  b 
*H 

u  «. 
o  « 
O'  o 


O'  4J 

C  -H 
•H  e 
U  M 
4)  C 
><  <0 
«  ^ 


I 


e 

c 

'0 


> 


a  O' 

fc  ^ 

o  < 


£ 

O' 


« 

u 


<0  c 

i3  O 


^  3 

a 

-»  o 


c  •-< 

«M  O  « 

O  O  ki 
^  *> 

•  ^  o  • 

>  10  V  u 

i3  a  c 
c  w  « 
E> 

O  £  T? 
«  C  O'  <0 
U  U  --4  u 
9  O  JZ 

•o  Xi  u 
c  u  xz 
3  -H  4>  -o 
O  9  to 


c  u 


c  > 

•H 

Ai  U 
><  £  < 
a  o* 
0*^0 
U  r1  « 

tf)  w 

o 


I 

.-3  O 
K  U 

z  *> 

•  o 
>«  o 
u  a 

O  M 
*J  •• 

2  E^'- 

O  O  I 
£  «-*  B 
<9  m  U 
J  C 


<0  o 


f>l 

41 

AJ 

c 

*0 

•H 

> 

CO 

la  B 

10 

>  C 

tf 

£  E 

£ 

la 

t 

C 

c 

•o 

CV 

41 

pa 

tf 

a  £  u 

tf 

« 

c 

« 

4J 

to 

40 

y 

£ 

10 

la 

pa 

(6  O 

(0 

3  «■« 

•pH 

•pH 

o 

GO 

> 

c 

a 

•»a 

C  •• 

la 

tf 

a  £  c 

£  O 

TJ 

la 

•pH 

C> 

1 

40 

o 

« 

la 

V 

*-4  » 

aJ 

£ 

-H 

£  tf 

0 

P)  • 

tf 

3 

£ 

> 

u 

« 

o 

0 

4) 

CO 

U 

tf 

£ 

O' 

w  o 

«p^ 

CC 

0 

3 

£ 

0 

u 

*0 

AJ 

M 

O  <^ 

41 

e 

p^ 

•H  £ 

3 

w 

« 

E. 

pH 

1 

1 

u 

«  O' 

c 

3 

U) 

O 

a 

o 

3 

^  B 

U 

£ 

U 

pa 

0 

4) 

o 

£  C 

•0 

£ 

I 

C  t 

V) 

la 

E 

4a  -H 

tf 

tf  C 

c 

tf 

T3 

m 

C 

*j  -a 

£ 

41  « 

tf 

tf 

ra 

pa 

-i  0 

« 

U 

tf 

« 

E 

O'  n 

E 

TJ 

£ 

O'  CO 

£  U 

•m 

0 

tf  «a 

£ 

in 

la 

* 

C  3 

£ 

•0 

•-I  (ft 

•  •• 

O  Cl 

O'  c 

U 

'« 

«  .. 

s 

U  £ 

£ 

£ 

tf 

H 

*0 

41 

AJ 

u 

<  c 

AJ 

>a  .H 

IQ 

V 

CO 

<-l  p-., 

tf  3 

•pa 

a  £ 

£ 

M 

c 

^  u 

-<H 

0 

M 

4J 

la 

£ 

£ 

CO 

£  «H 

tf 

a  ^ 

E 

o 

O* 

<0 

41  C 

la 

£  V) 

—a 

U 

U 

■0 

C 

<^ 

£  1 

u 

U  0 

in 

•o 

£ 

>  •-• 

0 

41 

O'  -a 

41 

*j 

..X. 

>  *0 

M 

•H 

3  E 

tf 

in 

c 

tf 

c 

£ 

m 

£ 

O' 

o 

•H  Wi 

c 

U  »J 

la  la 

t 

£  0 

la 

41 

tf 

^a 

3 

E 

ia 

> 

c 

S  <0 

V 

-a  o 

41 

c 

in 

U)  m 

H 

•«  la 

la 

.«a 

O 

>» 

£ 

» 

-•-< 

41 

40 

<0 

a 

40 

40 

E 

la  ba 

>  la 

o 

CO 

■ 

in 

£ 

£ 

la 

la 

T3 

e  a 

U  E 

—a 

»j 

3 

41  < 

••a 

(^ 

T3  O 

U 

CO 

— » 

O 

O' 

« 

T5 

<0 

E  « 

O'  £ 

•H  o 

A^ 

la 

JZ 

£  *0 

£ 

pa 

C 

•pa 

O'  tf 

S 

> 

ia 

O 

c 

£ 

C 

••a 

O'  u 

la 

A^ 

a 

A>  to 

3 

tf 

la 

pa  la 

C 

0 

Va 

T3  H 

-■H 

E 

41 

<0 

O* 

V> 

in  4J 

-a  C 

pa 

px. 

c 

tf 

£ 

0 

*0 

•• 

£ 

a  £ 
e  E 


u 

H 


K  -•  k- 


lilt 


U  £ 

C  O'  £ 

a  -*  e 


I  O  u 
£  U  U  O 
O'  W)  u  w 


£  --  U  • 

a  <  4)  < 
VI  (/)  a 
o  <  w  • 
E  Z 


U  f'  £ 


I  ia  V) 


^  C'  o 
-•  E 

.J 


/)?:a^urohCfit 


ta 

o 

a: 

o« 


o 

M 

n 

c 

a 

o 

o 

La 

c 

AJ 

4A 

u 

9 

*l>^ 

s 

to 

o 

VJ 

jC 

0) 

1 

4> 

4^ 

0 

3 

X 

0 

•-A 

«Q 

o 

3 

0 

•0 

0 

u 

c 

2 

tn 

<u 

a 

S 

*0 

a 

> 

3 

0 

0 

M 

o 

a 

rA 

X 

rH 

o 

i>* 

0 

o 

xz 


o  c 


4J  —* 

0  <u 
Q)  t: 
<u  o 

«u  E 
<b 

0  o 
£  c 

4J  IT] 


CT  TJ 

c  <a 

t* 

AJ  '•'v 
CJ 

u  o 
o  c 

CL  «a 

Lf  4J 

O  AJ 
0  -H 

c  m 

M  « 

c 

«3 
Li 


C  .H  >, 
— ^  M  ij 

I  -H  •  *0 
>«  Li  «  41 
i3  O  •  TJ 


O  C 
C  *»4 


> 

•'  o 
0)  u 
<11  a 

VI 

«  c 
cr  o 

•w 
e  4J 
«  3 
W  -H 
3  O 
0  (fl 
£ 

C  « 

«  -H 
«>  A> 
Li  Ll 
O'  « 

s  a 


X' 


0  O  0  o 


a. 

H 

a 

N 


n 

o 

i3 

«0 

-H 

W 

<0 

> 

4J 

« 

CL 

c*  m 

u  c 

Lt  « 

«  fl 

^  4) 

c  u 
3 
ei 
<0 
•  o 
^  a 


4i 

m 


cn 

Q 


«  a  c 

CU  JC  o 

<0  -P  •H 
O  -H 
^  Li  3 

c  o  ^ 
o  o»  o 

•H  w 

L>  « 
iQ  c 

iM  e  0 

c  a 

o  o 

-H  4J 

a  «  ^ 

U  La  a 
O  4] 
*  >  La 

TS  C  I 
La  -w  w 
«  K  t» 
? 

U  «  3 

o  >  cr 

4a  <«-#  C 

K  4> 

«  « 
*•  La 
9  0 
L»  4>  *0 
C  -H  -H 
«  > 
BOO 
«  4}  La 
La  La  0* 
3 

O  3  L> 

«  o*  o 

4)  4}  3 

a  La  a 


o  o  o 


c 

C7» 

•H 

0) 

4) 

■o 

4J 

c 

9 


La 

JJ 

M 

C 


n 

4; 

o 

>» 

w 

« 

c 

<Q 

4> 

> 

L> 

o 

c 

4) 

to 

O' 

c 


4> 

LA 

4) 

O 


O' 

c 


a 

3 

O 

o 

V 

c 


o 

4) 

LA 

4) 

a 

n 

/z 

to 

o 

c 


O' 

c 

•H 

•o 

c 

A 

L» 

o 

La 

0 

•o 

c 

:3 


c 

Q> 

AA 

o  c 
a  o 

•*A 

4]  a 
o  o 
c  c 

<4  *j» 


>•  4J  «, 


La 

<0 

3 

O 

•«A 

4L 

Li 

♦a 

a 

La  tn 

o  c 
*H  o 

>  •H 

10  O' 
l:  4) 
4>  Va 

O'  31 
C  0 
•H  TJ 
»  C 

La  > 
«  s 

4a 


a  u 

V)  < 
C  0 

10  c 

La 

•U  ^ 

I 

c  c 

•H  O 

c 

c 

o  c 

•^A  ^ 
D 

0  a 

3  0 

c  ^ 

4)  A 

«*j  o 

AJ  La 
IQ  a 

^  0 
^  o» 
IQ  La 
a  10 
U  w 


0 

> 


RAD:  A  LINE-BY-LINE  RADIATIVE  EXCITATION  MODEL,  WITH  APPLICATION  TO  CO2 
P.P.  Wintersi^incr*,  R.H.  Picard*,  R.D.  Shaima*,  H.  Ncbcl**, 

A.J.  Paboojian*,  J.R.  Winick*.  and  R.A.  Joseph* 

*ARCON  Corporation,  260  Bear  Hill  Road,  Waltham,  MA  02154 
*Geophysics  Laboratory/OPE,  Hanscom  AFB,  MA  01731 
^Permanent  address:  Alfred  University,  Alfred,  NY  14802 

In  the  non-LTE  region  of  the  atmosphere,  the  local  kinetic  temperature  is  not  sufficient  to  determine 
molecular  vibrational-lcvel  populations  (or  vibrational  temperatures)  as  it  is  in  the  LTE  atmosphere,  and 
it  is  necessary  to  calculate  these  populations  explicidy  in  order  to  determine  die  radiadon  fluxes.  This 
requires  a  calculadon  of  the  excitation  and  loss  processes  for  the  molecular  vibrational  levels,  which,  in 
turn,  requires  a  radiative-transfer  calculation  to  determine  the  upwclling/downwcUing  fluxes  responsible 
for  the  radiative  excitation  process.  We  have  developed  the  RAD  line-by-line  radiative-transfer  code  to 
carry  out  this  calculation  in  the  atmosphere.  This  paper  will  first  discuss  the  basic  RAD  algorithm  and 
approach.  Then  we  will  show  the  results  of  calculations  carried  out  for  CO2  15  fj.m  and  4.3  /im  radiation 
under  nighttime,  daytime,  and  auroral  conditions.  Finally,  we  will  validate  the  model  by  comparing  the 
RAD  predictions  to  data  from  the  SPIRE  and  FWl  experiments. 
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LINE-BY-LINE  RADIATIVE  EXCITATION  ALGORITHM 
(RAD)  -  FEATURES 


-  ASSUMPTIONS 

0  PLANE  -  PARALLEL  GEOMETRY 
0  NON-OVERLAPPING  LINES 
0  COMPLETE  FREQUENCY  REDISTRIBUTION 
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-  FULL  L'NE-BY-LINE  CALCULATION 

-  FINE  UYERING 

-  ALTITUDE-DEPENDENT  EMISSION/ABSORPTION 
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LINE-BY-LINE  RADIATIVE  EXCITATION  ALGORITHM 
(RAD)  -  APPLICATIONS 
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RAD:  A  LINE-BY-LINE  RADIATIVE  EXCITATION  ALGORITHM 

OUTLINE 
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-  APPLICATIONS/COMPARISONS  WITH  DATA 


.44j/v’  -  LTi>  hr  rir.  ec^r<-. 


iJ9Ly  T=6fl  PoPOL^  T/ cry's  V-  f' L  O^S^, 

0  'TdhN^fS'CZ.  -  Y^L,  U  X 


=  5  r 


SoVJftCg  fUf/Onor/ 

5t  — 


■?fV 


3.  jv 


i)gftl£^.  -  gw^.r/ov  (Sr£hw^^rArr) 

a/  ^X~  _  — 

^  ’  K.  +  A 


-  t>^f/A«r/<>yV  p-  T^.I 

K<  b 

-i^.-  , 

J  =  hvea^c^s  oP  X 


y 

t 

f 

■f 

1 

h  ®^| 

n 

/ 

1 

> 

1/  I 

j 

L. 

C.(?*-i-/$/W^  Rao/a 


RAD  MODULAR  STRUC^TURP 


LOCAL  AND 

RADIATIVE 

EXCITATION 


Populations 
or 

Vibrational 
Temperatures 


_Fluxes- 


RADIATIVE 

TRANSFER 


ITERATIVE  APPROACH  ->  SMALL  EQUATION  SET 

-  Method  of  Successive  Substitution 

-  Provisional  Vibrational  Temperatures 

— >■  Radiative  Excitation  Rates 
->  New  Vibrational  Temperatures 

-  Major  &  Minor  Iterations 

-  Minor  Iterations:  Opacities  Constant 
— >  Transfer  Function 


CO. 


2"  Co^  + 


B  J 

-  CouLi  $lOA/PrL-  (^'''^) 

Cojf+Al  ^  -^Al 

Xt 

-  i^i6ftftT)0AfAL  TfinMsrea  LV-v] 


1 0^  + 


Kty., 


+ AJ 


WV 


y'  ^  y 

^  3 


y  ^  y 
^  3 


CO2  BEND-STRETCH  MANIFOLD 


I2MI 


12M2 


M0Q1 

70007  "'  -■'■  "  I  0*401 

70003 


fIlOf 


ttt07 


03301 


6t7t.r,r‘ 


o  c.  m'  ■ 


I  I  !  '■ 


TANGENT  HEIGHT  (km) 


SPIRE  BAND  RADIANCE,  1 3-1  6.5 


RADIANCE  (watt/ cm^— ster) 


ALTITUDE  (km) 


VIBRATIONAL  TEMPERATURE  (K) 


I  - 


tangent  height  (km)  * 


200  ^  - -  j 

VIBRATIOKAL  TBIPERAWRE  W 


daytime  CO2  4.3  /zm  SPIRE 


lu  '  i(r 

RADIANCE  (W/cm2sr) 


C02  1^3 

(NIGHT  TIME) 


V1BRAT1(^4AL  TEMPERATURE  (K) 


CO2 1/3 

(NIGHT  TIME) 


180  200  220  240  260  280  300 

VIBRATIONAL  TEMPERATURE  (K) 


li'.MPtKATURE  (K) 


RAD  -  LINE-BY-LINE  RADIATIVE  EXCITATION  ALGORlTHf^> 


CONCLUSIONS 


POWERFUL  LINE-BY-LINE  RADIATIVE  EXCITATION  MODEL 
(RAD)  DEVELOPED 

RAD  APPLIED  TO  CO^  4.3^01  &  IS^im  AIRGLOW  PREDICTION 

MODEL  PREDICTIONS  AGREE  WELL  WITH  SPIRE  AND  FWI 
DATA  UNDER  BOTH  NIGHT  AND  DAY  CONDITIONS 

NON-LTE  EFFECTS  ABOVE  60-70  KM  (+/-) 

RAD  MODEL  PROVIDES  BENCHMARK  TO  TEST  VALIDITY  OF 
MODELS  WITH  MORE  APPROXIMATIONS  (HAIRM,  SHARC) 

NEXT  TWO  TALKS:  APPLICATION  TO  FASCOD3 

APPLICATION  TO  4.8^m  CO 
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NON-LTE  CO2  VIBRATIONAL  TEMTERATLRE  PROFILES  FOR  FASCOD3 
P.P.  Winicrslcincr  and  A.J.  Paboojian 
ARGON  Corp.,  260  Bear  Hill  Road.  Waltham,  MA  02154 
J.R.  Winick  and  R.H.  Picard 
Geophysics  Laboratory/OPE,  Hafiscom  AFB,  MA  01731 

Using  line-by-line  radiative  transpon,  we  have  developed  procedures  for  evaluating  non-LTE  popula¬ 
tions  for  all  the  infrared-emitting  vibrational  states  of  CO2  that  arc  important  for  atmospheric  studies.  Wc 
have  calculated  these  populations,  c.sprcsscd  as  vibraiiona!  temperatures,  for  tJie  six  FASCODE  model 
atmospheres  for  day  and  night  conditions.  Wc  will  outline  the  calculations,  including  the  assumption  and 
the  data  required  for  input,  and  present  some  of  the  results  (c.i*..,  vibrational  temperatures,  cooling  rates, 
.'xcitation/de -excitation  rates).  We  will  also  di.scuss  the  qualitative  differences  arising  from  the  different 
features  of  the  model  atmospheres. 


jn-LTE  CO.  Vibrational  Temperature  Profiles  for  FASCOD3 
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SUMMARY 


ALL  ABOUT  NON-L TE  CO2  IN  THE  ATMOSPHERE 

(or) 

WHERE  IS  NON-LTE  BEHAVIOR  SIGNIFICANT  ? 


IMPROVED  CO2  VIBRATIONAL-TEMPERATURE  LIBRARIES 

o  ALL  IMPORTANT  EMITTING  STATES 
o  UP-TO-DATE  MODEL 
o  LINE-BY-LINE  RADIATIVE  TRANSPORT 
o  CONSISTENT  WI  FI  I  FASCODE  ATMOSPHERES 


RESULTS 

o  SUBSTANTIAL  DIITERENCES  BETWEEN  MODEL  ATMOSPHFf 
o  DIURNAL  DEPENDENCE 
u  DEPENDF  NCE  ON  (Cn  .j,  (OI 
o  ISOTOPIC  DEPENIM  Ni:i 


NON-LTE  CO  INFRARED  EMISSION  IN  THE  MESOSPHERE  AND 
LOWER  THERMOSPHERE  AND  ITS  EFFECT  ON  REMOTE  SENSING 


J.R.  Winick,  R.H.  Picard 

Geophysics  Laboratory/OPE,  Hanscom  AFB,  MA  01731 

P.P.  Wintersteiner,  A.J.  Paboojian 
Arcon  Corporation,  260  Bear  Hill  Road,  Waltham,  MA  02154 

Wc  will  present  a  model  of  the  non-LTE  infrared  emission  from  CO  based  upon  the  RAD  computer 
code.  This  model  uses  a  line-by-line  radiative  traasfer  algorithm  to  calculate  the  radiative  contribution  to 
the  vibrational  temperature.  The  CO  emission  departs  from  LTE  at  altitudes  as  low  as  40  km  at  night  and 
30  km  during  daytime.  The  nighttime  CO  vibrational  temperature  above  60  km  is  most  strongly  dependent 
upon  the  radiative  excitation  emanating  from  the  upper  stratosphere  and  weakly  dependent  upon  the  CO 
density  profile  above  40  km.  This  introduces  an  additional  seasonal  and  latitudinal  dependence  to  the 
niglittime  radiance  beyond  that  caused  by  local  variations  in  the  CO  density.  The  daytime  vibrational 
temperature  is  dominated  by  direct  absorption  of  solar  radiation  which  has  important  structure  imposed 
on  it  from  solar  CO  absorption  lines. 
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CO  Non-LTE  Emission 


OUTLINE 


Review  of  GL  CO  infrared  Observations  (FWI) 

Use  of  line-by-line  non-LTE  model  (RAD)  for  determination  of  CO 
vibrational  temperature 

Important  Processes  Determining  Nighttime  non-LTE  CO(v=1) 

Sensitivity  Study  of  vibrational  temperature 

Comparison  of  CO  Vibrational  Temperature  to  Kinetic 
Temperature 

Important  Processes  Determining  Daytime  non-LTE  CO(v=1) 
Solar  Excitation 

Solar  CO  absorption  line  struciure  in  solar  flux 
Can  one  develop  retrieval  algorithms? 

Complex  fonvard  problem 


Conclusions 


KINETIC  MECHANISM  FOR  CO(V)  NLTE  EMISSIONS 

CO(v=1)  +  M  <==>  CO(v=0)  +  M  M=N2,02,  AE=2143  cm  ^ 

k^,=6.67x10'®exp(-208n'^'^) 

C0(v=1)  +  0  <==>  CO(v=0)  +  O 

k2^=1.56x10-''^exp{-118/T’^^) 

V-V  processes; 

CO{v=0)  +  N2(v=1)  <==>  CO(v=1)  +  N2(v=0) 
k3^=6.98x10*^^exp(-25.6/T^'^) 
k3, =8.42x1 0'^^exp(-48.5/T^'^) 

CO(v=0)  +  COjlOOl)  <==>  CO(v=1)  +  C02(000)  AE=20b  cm  • 
k^,=1.56ex10'^^exp(-30.in'^'^) 

CO(v=1)  +  O2(v=0)  <==>  CO(v=0)  +  02(v=1)  AE=636  cm^ 

k5j=3.5x10''°exp(-124/T^'^) 


AE=2143  cm’"' 


AE=188  cm'"' 


KINETIC  MECHANISM  FOR  CO{V)  NLTE  EMISSlON(Continued) 


Radiative  Processes: 

CO(v=1)  ==>  CO(v=0)  +  hv  =  34  s"* 

CO(v=0)  +  hv  ==>  CO(v=1)  Jco 

Calculation  of  Excitation  coefficient  in  NLTE  atmosphere  done 
by  line-by-line  Code  RAD 

In  the  mesosphere: 

V-T  is  negligible  <  10'^  s'  ,  k^^[M]  «10  s 

V-V  with  O2  (Reaction  5)  is  unimportant 
kg  is  dominant  collisional  term  and  depends  upon  N2 
vibrational  temperature 

Reactions  4  is  unimportant  since  CO^  is  a  minor  species 
Reactions  with  atomic  oxygen  only  important  at  high 
temperatures,  but  above  150  km  density  Is  too  low 


RADIATIVE  TRANSPORT  -  KINETICS  CALCULATION 

Line-by-LInc  Plane-Parallel  Atmosphere  Code 

Lineshape  changos  v/iih  altitude:  Voigt  lineshape 
Doppler  broadening  (temperature) 

Lorentz  broadening  (pressure) 

AFGL  HITRAN  Database  using  79  (1-0)  lines 

Vibrational  temperature  calculated  by  iterative  method 
Steady-state  calculation  of  radiating  state 
Uses  previous  iteration  radiation  field 
Combines  collisional  excitation  and  quenching 
CO  profile  is  unusual  in  that; 

Mixing  ratio  increases  above  Tropopause 
Latitudinal  variation 
Seasonal  variation 
Collisional  processes  are  inefficient 


CO  PROFILES 
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CO  Vibrational  Temperature:  SZA  Dependence 
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CO  Non-L  TE  Emission 


Conclusions 


Nighttime: 

CO  4.8  |j,m  nighttime  Emission  is  non-LTE above  45  km 

Collisions  are  inefficient  populating  CO(v)  which  leads  to  the  dominance 
of  radiative  processes  above  50  km 

CO  has  unusual  profiles  that  make  it  a  good  test  for  radiative  transfer 
codes 

Latitudinal  and  seasonal  variation  in  CO  and  T  effect  ^^{[3(2) 

Large  departures  from  LTE  can  lead  to  large  errors  in  retrieved  density  if 
non-LTE  model  is  not  used 

Limb  radiance  retrieval  is  especially  difficult  since  the  dependence  upon 
the  lower  state  density  is  partially  non-local 

This  non-local  problem  may  not  be  too  important  except  for  mid 
and  high  latitude  winter  situations  where  there  is  increased 
CO  in  the  lower  mesosphere 


CO  Non-LTE  Emission 

Conclusions  (continued) 


Daytime: 

CO  4.8  pm  daytime  Emission  is  non-LTE  above  30  km 

Production  of  CO(v=1)  is  dominated  by  direct  solar  excitation 

CO  solar  absorption  lines  must  be  included  in  excitation  calculation 

Daytime  enhanced  N2{v=1)  and  CO2(001)  populations  effect  excitation 

(only  weakly,  rate  constant  and  SZA  dependent)  in  60-80  km  region 

Vibrational  Temperature  profile  depends  upon  SZA  especially  for  SZA>75 

Mesospheric  Emission  is  less  sensitive  to  latitude  and  season  (except  for 
S2IA  dependence)  than  nighttime  case. 

Subarctic  winter  CO  profile  has  less  solar  excitation  due  to 
slightly  increased  opacity  at  high  SZA. 
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PATH  CHARACTERIZATION  ALGORITHMS  FOR  FASCODE 


R.G.  Isaacs,  S.A.  Oough,  R.D.  Worsham,  J.-L.  Moncet,  B.L.  Lindner,  and  L.D.  Kaplan 
Atmospheric  and  Environmental  Research,  Inc.,  840  Memorial  Drive,  Cambridge,  MA  02139 

We  describe  the  results  of  a  study  undertaken  at  AER  to  identify  and  implement  a  state-of-the-art 
nonlinear  retrieval  approach  to  characterize  line  of  sight  variability  of  atmospheric  thermal  and  constituent 
environments.  This  path  characterization  capability  was  designed  to  interface  with  the  existing  Geophysics 
Laboratory  (GL)  line-by-line  radiancc/transmittance  code,  FASCODE. 


aer 


Path  Characterization  Algorithms  for  FASCODE 

R.  G.  Isaacs,  S.  A.  Clough,  R.  D.  Worsham,  J.-L. 
Moncet,  B.  L.  Lindner,  L.  D.  Kaplan 


Atmospheric  and  Environmental  Research  Inc. 
Cambridge,  MA  02139 


aer 

Path  Characterization 


Quantitative  assessment  of  the  relationship  between  path  optical 
properties  and  path  thermodynamic  properties; 

•  Path  optical  properties  (data): 

•  optical  depth 

-  transmitttance 

-  emission  (radiance,  brightness  temperature) 

•  Path  thermodynamic  properties  (desired  parameters): 
-temperature 

-  concentrations  of  absorbers  and  scatterers 

-  path  boundary  temperature  and  emissivity/reflectivity 

-  pressure  (surface,  limb  layer) 


Objectives 


•  Identify  and  implement  a  state-of-the-art  nonlinear  retrieval 

approach  to  characterize  line-of-sight  variability  of  atmospheric  path 
thermal  and  constituent  environments; 


•  Interface  with  the  existing  Geophysics  Laboratory  (GL)  line-by-line 
radiance/transmittance  code,  FASCODE 

•  Provide  code  and  comprehensive  documentation 


aen 


Path  Characterization  Concept 


Select 


Simulate 


Invert 


Diagnose 


Screening  Algorithm 


aen 


•  Suggest  fruitful  spectral  regions  for  the  retrieval  of  desired  atmospheric 
parameters; 


•  User  selects  temperature  parameter  or  molecule  of  interest; 


•  Algorithm  computes  optical  depth  vs  total  optical  depth  for  selected 
wavenumber  region; 

•  HITRAN  data  base,  FSCATN  atmospheres  (T,  p,  u  ),  1  cm'^ 
band  model; 


•  Histogram  capability  to  select  wavenumber  regions  for  best 
discrimination  from  background  atmosphere. 


aer 

Forward  Problem  (FASCODE) 


•  Provides  simulated  data  set  (radiance  spectrum,  channelized 
brightness  temperature  data,  etc.)  from  initial  guess; 


•  Provides  required  Jacobian  of  data  sensitivity  to  variations  in 
desired  parameters:  p  . 


dx  . 

J 

•  Covers  microwave  to  UV  spectral  domain; 


•  Interfaces  provide  all  required  FASCODE  data  files  in  standard 
format; 


•  Capability  to  control  data  spectral  resolution; 


Computational  liability  potentially  addressed  via  "rapid"  algorithms 
for  selected  channel  sets  . 


Retrieval  Algorithm  Selection 


•  Literature  review  (Isaacs,  1988,  SPIE,  928, 136)  prompted  choice  of 
Physical  Least  Squares  (PLS)  approach; 


•  Extensive  heritage  (Rodgers,  1976, 1987;  Eyre,  1989); 


•  Incorporates  the  physics  of  the  radiative  transfer  process  through 
the  FASCODE  forward  problem; 


•  Climatology  not  required  (helpful  In  practice,  null  space,  1st  guess) 


•  Applied  iteratively; 


aer 


Retrieval  Algorithm  Capabilities 


•  Choice  of  penalty  function; 

-  maximum  likelihood  (using  error  covariances  of  parameter  and  data) 

-  ridge  regression  (minimum  information); 

•  Choice  of  representation:  radiance,  brightness  temperature; 

•  Option  for  eigenanalysis  including  retrieval  via  eigenvalues; 

•  Capability  to  include  physical  constraints  e.g.  superadiabaticity/ 
supersaturation,  etc.  (not  implemented); 

•  Capability  to  add  realistic  measurement  noise  to  simulations. 


aer 


Error  Analysis 


•  Provides  comprehensive  error  analysis  based  on  Rodgers  (1987); 

•  Estimate  of  the  covariance  of  the  retrieval  errors  including  the 
effects  of: 

-  measurement  noise 

•  uncertainties  in  the  a  priori  information; 

•  Errors  in  forward  problem  not  yet  treated; 

•  Retrieval  performance  evaluated  by  comparing  covariances  of  the 
state  parameters  before  and  after  the  measurement  process: 

-  information  content 

-  FUV 
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Test  Cases 


•  Microwave  temperature  and  moisture  sounding: 
Advanced  Microwave  Sounding  Unit  (AMSU); 

•  Infrared  temperature  and  constituent  (ozone)  sounding: 
SCRIBE  and  HIS; 


*  Spectroscopic  applications. 


aer 

Path  Characterization  Summary 


•  implementation  of  a  state-of-the-art  nonlinear  retrieval  approach  to 
characterize  iine-of-sight  variability  of  atmospheric  path  thermal  and 
constituent  environments  consistent  with  the  FASCODE  formalism; 


•  technical  issues  reiated  to  non-linearities,  incorporation  of  physical 
constraints,  forward  problem; 


•  applications  to  other  path  parameters  such  as  aerosol,  cloud,  and 
non-LTE. 


•  Final  Report:  Isaacs  et  al,  "990  (GL-TR-90-0080) 


TEMPERATURE  RETRIEVALS  WITH  SIMULATED  SCRIBE  RADIANCES 


J.-L.  Moncct,  S.A.  Clough,  R.D.  Worsham,  R.G.  Isaacs  and  L.D.  Kaplan 
Atmospheric  and  Environmental  Research,  Inc.,  840  Memorial  Drive,  Cambridge,  MA  02139 

The  FASCODE  path  characterization  algorithm  has  been  applied  to  the  retrieval  of  atmospheric 
temperature  using  simulated  radiances  for  a  nadir  view  from  31.3  km  corresponding  to  a  typical  SCRIBE 
float  altitude.  The  spectral  region  studied  was  720  to  775  cm  ’  with  a  resolution  consistent  with  SCRIBE. 
The  sensitivity  of  the  retrievals  to  measurement  noise,  first  guess  error,  and  to  systematic  error  including 
photometric  calibration  and  line  parameter  error  will  be  discussed. 


Temperature  Retrievals  with 
Simulated  SCRIBE  Radiances 


J.-L.  Moncet,  S.A.  Clough,  R.D. 
Worsham,  R.G.  Isaacs  and  L.D.  Kaplan 


Introduction 

Goal  of  the  study; 

Assess  the  pcrfonnance  of  path  characterization  algorithm  for  temperature 
retrievals  based  on  SCRIBE  data. 

Study  effects  of; 

•  instrumental  noise 

•  calibration  errors 

•  uncertainties  on  line  parameters 


SCRIBE  Instrument 


•  Cryogenic  Michelson  interferometer  spectrometer 

•  Platfonn:  Balloon 

•  Altitude:  32km 

•  Spectral  range  600-1400  cm“^ 

•  Resolution:  0.06cm~^ 


Non-Linear  Retrieval  Algorithm 

Obtain  best  estimate  x  of  unknown  state  vector  from  set  of  measured 
radiances  y"\ 

y- =  F  (X-) -i- s  , 

where  e  (covariance  S^-)  is  the  measurement  error. 

Newtonian  Iteration: 

•  Expand  forward  model  about  guessed  vector  Xn-i- 

y  Vn—i  K(a:  Xn_i)  “bO(x  .  {2} 

•  Ignoring  higher  order  terms  in  eq.(2),  find  Xn  such  that 

E{xn)  =  -  yn)”^  W(y'"  -  y„)  -  Xo)’^r(x„  -  Xq)  (3) 

is  minimized  (here,  W  —  S"')- 

•  Solution*: 

Xn  =  atn-i  +  H'*  [K'^W(y^  -  y„_i)  +r(xo  -  a^r.-i)]  (4) 

with 

H-^  =  (K'^WK  Er)"'.  (5) 


*  In  Maximum  Likelihood  approach,  T  =  bj  ‘ ,  where  Si  is  the  covariance  of  the  first  guess. 


Error  Analysis  (Rodgers,  1990) 

•  Once  convergence  is  readied,  and  provided  that  K  does  not  change  too 
rapidly  within  the  error  bounds  of  the  solution, 

U"'  -  Vv.-i  —  K  (x’’‘  -  x,i_i)  +  e  .  (6) 


•  Combining  equations  (/))  and  (4),  the  retrieval  error  is  expressed  as  the 
sum  of  a  "null-space”  error  and  the  error  due  to  instrumental  noise; 

X..  -  -  H  '  nCx.,  -  x”‘j  eR-^K'^Ws  .  1- 


Ehe  total  error  covari.ince  is. 


with 


S  "  S  • '  -  s 

- 11  rsd'ir’ 
iid-ow'b.wKir 


and 


Conditions  of  the  Experiment 

•  Temperature  retrieval  in  the  C02  15/rm— band  (720  -  775  cm”') 

•  Model:  FASCOD3 

•  Line  tile:  HITRAN  1988 

Atmospheric  Profiles: 

•  Retrieved  parameters:  18  layer  temperatures  (from  0  to  30.311cm)  + 
surface  temperature 

•  Temperature  profile: 

-  Background:  Tropical 

-  Simulated:  Background  +  2K 

•  Water  vapor  profile:  U.S.  Standard 


ATMOSPttERIC  PROFILES 


- TROPICAL 

—  US.  STANDARD 


;  t 


'  f'  )  I 


'  I  I 


*Tc«vf4ro.'uf«.  84^ri«,v*.\*  S  cf»"  ) 

t,rror&  on  tin*.  P^rn.mo'Vtr  i 


TEMPEHATUBE  DIFFEHENCES  (K) 

♦  f  ft*«)Ua*«**  u>.‘4U 
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Summary 

•  Retrieval  most  effective  in  lowest  atmosplieric  layers  (below  5km)  and  at 
the  surface.  Above  5km  tlte  improvement  depends  on  the  quality  of  the 
a  priori  information  (if  short-range  forecast  is  used,  little  improvement 
expected  in  the  10-30km  layer). 

•  In  an  ideal  situation  (no  model  errors)  expected  accuracy  is  <  1-2K  in 
the  lowest  lOkm,  depending  on  level  of  instrumental  noise. 

•  Retrieval  is  very  sensitive  to  calibration  errors.  These  tend  to  affect  all 
layers,  incluoing  the  surface. 

•  Current  state  of  knowledge  on  line  paraineters  seriously  limits  the  quality 
of  the  physical  retrieval  (ne.ir  the  surface,  retrieved  temperatures  remains 
little  affected  bv  errors  on  line  intensitie.s). 


VALIDATION  OF  RANDOM  PHASE  SCREENS 
A.E.  Naiman  and  T.  Goldring 
W.J.  Schafer  Assoc.,  Inc.,  Arlington,  VA  22209 

A  simple  scheme  for  validating  realizations  of  atmospheric  turbulence  phase  screens  is  presented  and 
cntiqued.  The  stati.stics  of  interest  arc  the  first  few  Zcmike  modes.  For  screen  generators  based  of  the  Fast 
Fourier  Transform,  computed  and  theoretical  variances  ((!])  are  compared,  and  the  discrepancy  between 
them  is  analyzed.  A  method  for  quantifying  the  accuracy  of  Zcmike  mode  variances,  given  grid  size  and 
resolution,  is  discussed.  Finally,  numerical  results  are  presented  for  a  specific  turbulence  spcctmm. 

Reference; 

( 1 1  R.J.  Sasicla,  A  Unique  Approach  to  Electromagnetic  Wave  Propagation  in  Turbulence  and  the  Evalu¬ 
ation  of  Multiparameter  Integrals,  Lincoln  Laboratory,  Massachusetts  Institute  of  Technology,  Technical 
Report  807,  1  July  1988. 


Validation  of 
Random  Phase  Screens 


Annual  Review  Conference  on 
Atmospheric  Transmission  Models 
5-6  June  1990 

by  A.  E.  Naiman  and  T.  Goldring 
W.  J.  Schafer  Associates 
Arlington,  VA  22209 
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A  Simple  Validation  Scheme 


•  Create  ensemble  of  phase  screens 

•  Extract  desired  Zernike  modes  directly  from  phases 

•  Eistimate  variances 

•  Compare  with  vaiues  expected  from  theory 


Estimated  and  theoretical  Zernike  mode  variances 
can  be  compared  directly  —  with  inherent  limitations,  . . . 

ATM  Conf." 


Validation  —  Definition 


•  Another  validation  criterion: 

•  Comparing  with  structure  function 

•  Zernike'  mode's  ri'pre’sent  classical  ejpticai  al)erratiejns,  e.g.  tilt,  focus 

*  lest  if  a  [)h<ise  screen  generator  realizes  the:  Zernike  moeles  pro[)eri\ 


AIM 


Purpose  of  Presentation 


•  Discuss  inherent  limitations  of  simple  validation  scheme 

•  These  arise  due  to  discretization 

*  Discrepancies  between  numerical  and  theoretical  results 


•  Ana’yze  these  errors 

*  Show  method  for  quantifying  them 
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Zerilike  Expansions 


Zerilike  expansion  of  a  continuous  2-D  function  <p  [Noll]; 

4>{r,0)  =  ^ajZj{p,e) , 

} 

R  -  aperture  radius,  p  =  '  €  [U,  1],  and  the  cocflicients  are  given  by: 

^  II  p  dp  (id  w{p)<f){Rp,e)  Zj{p,e) , 

I  *  ^  r  1 , 

\V{r)  =  ^’  ~  ’is  the  radially  symmetric  step  function 

0,  r  >  1, 


Any  continuous  2-D  phase  realization, 
can  be  expanded  in  terms  of  Zernikc  polynomials 
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Variances  of  Theoretical  Zernike  Coefficients 


Zernike  coefficient  phase  /ariance  [Sasielaj: 

==  JJdfl  UV(«)  C(/?) , 

U'^(«)  is  the  turbulence  spectrum,  Fj{k)  is  the  filter  function: 

2cos^(m'(p)  (x-component), 

2siir(mV)  (y-com[)onent), 

1  (m'  n), 

1)  diaiiirler  of  the  aperture, 

rn\  ri'  Z('rnik(’  component  indices  (ang\ilar  frequency  and  radial  degree 


k) 


in'  -  1) 


2,/., 


*11  2  I 
2 


I'. Ik'  ,./( Zernike  {lolynoniial  j)| 


N.,:i 


Effective  vs.  Tiieoretical  Variances  — 
Definition 


•  cr^|f(a_|)  variance  of  Zernike  coefficients  of  phase  screens 

*■  i.e.,  not  of  realizations  of  the  continuous  random  surface,  which  exist 
only  in  theory 

•  ^  standard  statistical  estimate  from  a  finite  sample  of 
phase  screens,  => 

*  cTesti^j)  ^  sample  size  increases 

•  However,  a^j(aj)  /  ^  fixed  computational  grid 

•  For  FFT-based  methods,  we  quantify  this  error 
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FFT  Phase  Screen  Generator 


•  For  grid  leugth  L,  with  N  points  in  each  direction,  our  phase  realization  is 
given  by  an  inverse  Discrete  Fourier  'IVansform  (DFT): 

0(l.!/)  =  EZ  6, 

n.Tii—  -  A 

a  continuous,  doubly-periodic  function  of  x  and  y. 


•  The  inverse  FFT  efficiently  computes  the  DFT  on  the  grid 
■  i,A:  =  0, ..  .,iV  -  1}  where: 


-  Vj  = 


JL 

N' 
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FFT  Phase  Screen  Generator  —  Continued 


•  'the  Fourier  coefficients  6„,„: 

‘  Random  zero-mean  Gaussian  complex  numbers 

•  <  h'cur  in  conjugate  p.iirs  -> 

+  (p(-c,?/)  is  real 

*  1  lav('  variance  of; 


Note 


'ncinsion  oi  a.v 


A  I  M  <  '•  aiferein  r 
p  i: 


FFl’  Geiier;itor  Choice  of 


a;o  determines  periodicity  and  frequency  content 

Choose  cJo  —  It  so  that: 

+  the  period  is  2L 

*  Extracting  any  contiguous  N  x  N  submatrix  of  y^)}  => 

*  Pairwise  covariances  will  be  correct  for  separation  distances  <  L 

*  The  continuous  random  surface  <^(x,  y)  is  sampled  at  precisely  the 
Nyquist  frequency  corresponding  to  the  bcind  limiting  assumption 

^  0,  j/t]  >  Nuj{). 

This  is  implicit  in  the  discretization  procedure. 


These  realizations  (p(x,  y)  are  distinct  from  the  theoretically  given  ensemble 
4>{x,y)  —  fj  whose  members  are  not  periodic. 

-'X'. 
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Effective  Zeriiike  Coefficients 


VVi'  have,  in  general,  the  Zernikc  coefncienls: 

Qj  -  II  pdpdO  \V{p)(p{lip,0)  Zj{p,0) , 
and  for  the  Ft  T  method,  the  phasi'  screen: 


rt  .T/t  :  :  -  A' 


or  e(}uival('utlv:  (p{r)  v  6,„„e' 

n  .m  •  -  A 


where  r  (  r.  jy)  and  --  ( 7iw,i,  ma’u).  Substituting: 


s  -  vy  h,,„jidd\V(p)Zj{p,o)c''^-''^ 

u  ui  A 


1  he  iiiiejqal  is  the  I'Miirier  t  ransform  of  the  Zernike  polynomial, 
over  iht'  aperture,  (ompuleil  analytic  allv  |Nollj  as  the'  tiller  function. 


A  1  M  (k>nfrr*‘fi  • 


mgn^ 

Tlieoretical  vs.  Effective  Variances 


Computation 


•  Taking  the  variance; 

n,7n—-N 

N-l  .y 

—  EE  C(^fim) 

n,m=-N 


•  Conclude:  Effective  variance  is  a  standard  numerical  approximator  to 
theoretical  variance,  i.e.: 

<Tl(a,)  =  lldKWt{K)F,(K), 

jV-1 

~  EE  4)(^^nin)  ^elT(^j) 

n,m=  — ,V 


is  “rectangle  rule”  approximation  to  cr^faj) 
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Theoretical  vs.  Effective 
Variances  —  Continued 


•  In  one  dimension: 


•  Note  the  error  d('pendeurie.s; 

‘  "rectangle"  error  |  as  a;,)  |  (  L  i) 
•  tail  error  [  <is  .V  I 


If  u.’,',’  is  not  included  in 
tin'll  (7,':,j-(uj)  *  0(.:,  as  /,,j\'  •  oc 


•V  I  M  l',.ul.-frn.  ' 
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Example:  Tilt  Variances 


•  In  general,  two-axis  tilt  variance; 


due  to  averaging  the  tilt  phase  variances  over  the  aperture 


•  Theoretical  and  effective  tilt  variances: 


and 
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Turbulence  Spectrum 


2-D  von-Karman  phiise  spectrum  [Tatarski]; 

^  0.033 

=  27r  fc' A.  -  , 

[4  +  /c2]  6 


where: 


k  =  the  spatial  wavenumber, 

A 

Ar  =  the  slab  length, 

cl  =  the  structure  constant, 

27r 

«o  =  ,  and 

Lq  —  the  outer  scale  of  turbulence. 


ATM  CoaStTtiice 
p  19 


Comparison  of  Theoretical  and 
Effective  Tilt  Variances 
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Error  Between  Theoretical  and 
Effective  Tilt  Variances 


•  Introduction 

•  Theoretical  and  Estimated  Zernike  Modes 

•  Error  Estimation  and  Effective  Zernike  modes 

•  Numerical  Example 

>  Error  Analysis  Theory 

•  Summary 
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Error  Analysis  for  2-D  Rectangle  Rule 


=  W4>{k)F{k)  ^ 

Rim  -  |«  :  |n  -  cl;o  <  +  -j  ‘^o.  -  - j  <  «:y  <  +  2)  ^0} 

<3  =  js:  -(w+2)‘^o<k.,k,  <(N-i)}=__^U_^P„m 

•  Expand  g(Kz,Ky)  in  a  Taylor  Series  about  the  center  point  {nuJo,mujQ)  of 
P„„,  and  integrate  term  by  term;  valid  for  n,  m  not  both  zero 
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Integrand  for  Small  Values  of  k 


In  (l  +■  I0‘'g( 


Small  Oscillations  in  Integrand 


l)oth  //  and  V7  arr  functions  of  k  r-  |/t]  alone 


Resultant  Error  Term 


Now,  5(0,0)  =  F(0)  =  0,  => 

N  —  \ 

^eff(aj)  =  ujlg{nujQ,mu}Q) 

ti  m  =  -  /V 
/  0 

Therefore,  neglecting  6^  and  higher  order  terms: 

aJCoj)  -  (T^ff(aj)  «  //  g  (Lk^  dKy  Jj  g  d/Cj  dKy  + 

Po.o  Q 

CJq  /V-1 

-  E  Y.  Vg{nujQ,  muja) 

n,  m  =  -A' 
n^+  m’  /  0 

g  is  circulajly  symmetric  => 

*  Double  integrals  can  be  approximated  by  single  integrals  in 
K  =  ^kI  +  kI 
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Relative  Magnitude  of  Error  Terms 
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Summary 


•  Meth(jd  for  validating  phase  generators; 

»  C'oniputf'  Zernike  modes  from  physical  phase  values 
»  t’oni[iare  with  continuous  and  discrete  theories 

•  nisrrepancv  between  two  theories  is  intrinsic  to  computational  necessity  of 

using  ,i  dis<  rete  grid 

•  Results  computed  from  actual  [)hase  screen,  at  best  rnatcli  up  with  discrete 
iht'cirv 
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Suiniiiarv  Coiitiiiuod 


•  Method  for  quantifying  accuracy  of  Zernike  mode  variances,  given  grid  size 
and  resolution 

•  Numerical  results  for  the  case  of: 

*  FFT  phase  screen  generator 

*  von-Karman  spectrum 

*  Tilt  variance 
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RECENT  ADVANCES  IN  MODELINC;  OF 
BOUNDARY  LAYER  REFRACTIVITY  TURBULENCE 

V.  niicrmann  and  A.  Kohnlc 
FGAN-Forschungsinsiitul  fur  Optik 
Schloss  Krcssback,  D-7400  Tubingen,  F.R.  Germany 

Monin-Oboukhov  similarity  is  an  elegant  tool  for  estimating  structure  parameters  and  inner  scales 
of  turbulent  refractive  index  fluctuations  in  the  boundary  layer.  Our  recent  experiments  support  slightly 
different  scaling  expressions  than  commonly  used.  Thc.se  modified  expressions  arc  in  better  agreement 
with  tlicorctical  considerations  and  lead  to  up  to  30%  larger  structure  parameters  and  up  to  15%  larger 
inner  scales. 

Tfic  program  SPIRIT  is  prcsci'.ied.  It  computes  structure  parameters  and  inner  .scales  in  tlic  lower 
boundary  layer  over  land  from  oasy-to-c.siiinatc  input  parameters.  A  listing  of  the  required  input  quantities 
and  an  example  of  model  results  in  comparison  with  measurements  arc  given. 
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rar;niieters  describing  the  spatial  spectrum  of  refractive  index 
in  the  optically  relevant  range; 

spectral  amplitude  in  the  inertial  subrange  (mm's  to  m's) 

-•  important  for  ail  turbulence  sensitive  EO-systems 

high  spatial  wavenumber  cut-off  of  inertial  subrange  (mm's) 
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computes 

Structure  I’arameters  and  Inner  scales  of  Refractive  Index  'I'ur bulence 
for  optical  frequencies  in  tiie  lower  boundary  layer  over  land. 


.Structure  of  SPlRI  f; 

1.  Assistance  for  estimating  input  parameters  (optional) 

2.  rarametciication  of  turbulent  fluxes  of  beat,  moisture 
and  momentum 

3  A!>!)licatiun  of  Monin-Oboukhov-Similarity 

■I.  I'tclatinj;  temperature  and  Immidity  to  refractivity 


INi’i;  i'  PARAMf- ri;KS 


Assistance  niodule; 


-  date,  time 

-  geographical  position 

-  height  above  sea  level 

-  cloud  cover 

-  wind  velocity  at  a  certain  height 

-  coverage  of  ground  by  obstacles  and  vegetation 

-  height  of  obstacles  and  vegetation 

-  humidity  of  the  ground  (humid,  moderately  humid,  dry,  very  dry) 

-  brightness  of  the  ground  (4  options) 

-  air  temperature 

-  (relativ  humidity) 

-  wavelength 


direct  entering; 

-  solar  irradiance  (day) 

-  strong  wind  heat  flux  (night) 

-  surface  roughness  length 

-  wind  velocity  at  a  certain  height 

-  ground  humidity  parameter  (0...1) 

-  ground  albedo 

-  air  temperature 

-  temperature  dcrivn.jve  of  refractivity 

-  humidity  derivative  of  refractivity 


MODfIL  OUTPUT 
Pnd'ile  of: 

-  struct  Life  parameters  C)". 
■  mn-r  scales  r,, 


mm  heinjit  de{)endr,  on  conditioms,  day;  dOO  iii 


S  y  ■>  : 


<  cA  r  >/o  ;  t/vA.  <*  t/f  ( 


'fJ'V.t.  (u(S7j 
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VISIBILITY  DATA  FILTERS  FOR  EUROPE 

B.A.  Schichlcl  and  R.B.  Husar 

Center  for  Air  Pollution  Impact  and  Trcitd  Analysis,  Washington  University,  St.  Louis,  MO  63130 

llie  purpose  of  the  report  is  to  present  the  mclliodology  for  filtering  die  European  meteorological 
visibility  data  from  undesirable  and  erroneous  data.  Seven  data  filters  were  devised  and  impo.scd  on 
the  European  synoptic  visibility  data  set.  The  data  .set  consisted  of  fourteen  years  of  meteorological 
data  (1973-1986)  for  about  16(X)  stations  in  Europe.  The  European  data  set  was  extracted  from  the 
DATSAV  global  weather  database  maintained  by  the  U.S.  Air  Force,  ETAC,  Sailt  Air  Force  Base.  31ie 
raw  meteorological  data  .set  consisted  of  over  1000  magnetic  tapes  containing  about  30  gigabytes  of  data. 
The  first  step  in  die  data  processing  involved  compacting  the  data  .set  into  a  binary  fonn,  which  reduced 
the  data  si/c  to  3  gigabytes.  Next,  from  the  daily  visibility  data,  the  quarterly  cumulative  distribution 
functions  for  extinction  coefficient  B,.,,  and  vi.sibility  were  computed.  Most  of  tlic  subsequent  data  filtering 
was  perfomicd  using  the  distnbution  functions  and  the  Voyager  data  exploration  software.  The  resulting 
database  is  suitable  for  in[)ut  to  radiative  transmi.ssion  and  transfer  models,  global  climate  models,  air 
pollution  studies,  as  well  tis  to  global  biogeochemical  explorations. 
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Figure  17.  Localioa  of  llie  1600  tUtioai  before  (he  filleri  were  upplicj. 
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Figure  6.  Coolour  maps  ol  cxlinclion  coefficient  (10*^  “*')  unlillcrcd  daU 
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COMPARISON  OF  XSCALE  89  VERTICAL  STRUCTURE 
ALGORITHM  WITH  FIELD  MEASUREMENTS 


R.P.  Ficgcl 

U.S.  Army  Laboratory  Command,  Atmospheric  Sciences  Laboratory, 

White  Sands  Missile  Range,  NM  88002 

The  visible  extinction  versus  altitude  data  from  the  Meppen  1980  and  Cardington  1984  balloon 
flight  field  tests  are  compared  to  the  recent  modifications  in  the  below  cloud  (Case  2)  vertical  profile  of 
XSCALE  89.  The  visibility  and  humidity  data  from  the  Sprakensehl  1983-1985  tower  measurements  are 
also  compared  to  the  low  visibility  (Case  1),  below  cloud  (Case  2),  and  radiation  fog  (Case  3)  model 
profiles  of  XSCALE  89.  While  certain  days  are  not  well  modeled,  we  find  good  agreement  overall. 
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the  low  cloud/fog  t:iscs.  The  line  has  slope  1.  top  is  at  7G  na. 


the  chvtuJ  top  is  at 


nas  Slope  l,  sligiiUy  more  ;  .ts  arc  within 
error  bar  of  agreement  tha  are  farther  away. 
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tlir  (lata  is  70  ±  7  from  the  model  73  %. 


SPATIAL  NON-UNIFORMITY  OF  AEROSOLS  AT  THE  15^00  ft  LEVEL 


L.A.  Mathews 

Naval  Weapons  Center.  Code  3892,  China  Lake,  CA  93555 
P.L.  Walker 

Physics  Depaitment,  Naval  Postgraduate  School,  Monterey,  CA  93943 

Longjump  III  meteorology  was  reported  on  previously*.  This  work  presents  the  results  of  a  more 
detailed  analysis  of  the  aerosol  data.  That  data  was  obtained  from  aircraft  at  15,500  feet  over  a  100  mile 
course  between  mountain  peaks.  Samples  were  taken  at  three  minute  intervals.  In  general  more  haze  was 
associated  with  the  mountain  peaks  than  for  the  region  between  them.  Extinction  fluctuated  at  shorter 
intervals  (10  miles)  along  the  flight  path  possibly  due  to  the  uplifting  action  of  convection  cells,  but 
superimposed  on  this  is  noise  caused  by  the  short  sampling  intervals.  Visibility  was  highly  nonuniform 
being  dominated  by  molecular  scattering  punctuated  by  regions  of  lower  visibility. 

‘Characterization  of  the  Atmosphere  for  Longjump  HI.  Mathews  and  Walker,  IRIS,  Ames  Research 
Center,  16  March  1989. 


SPATIAL  NON-UNIFORMITY  OF  AEROSOLS 
AT  THE  15,500  FOOT  LEVEL 


by 

l.A.  Mathews 
Research  Department 
Naval  Weapons  Center 
China  Lake,  Calif.  93555 

and 

P.L.  Walker 
Physics  Department 
Naval  Postgraduate  School 
Monterey,  Calif.  93943 


PROJECT  LONG  JUMP  III 


•  FIELD  TEST  OF  VARIOUS  DEVELOPMENTAL 

INFRARED  INSTRUMENTATION 

•  CONTROLLED  AIRCRAFF  TARGETS 

•  SPECIFIC  FLIGHT  PROFILES 

•  1 1,500  TO  15,500  FP 

BARCROFT  LABORATORY  (ELEV.  12,470  FT) 

•  ABOVE  MAJOR  PORTION  OF  ATMOSPHERIC  CONTAMINATES 

•  PROVIDES  A  HORIZONTAL,  CL.EAR  LINE  OF  SIGHT  FOR 

OVER  100  MILES 


TO  DETERM INfl; 

(1)  OPTICAL  PROPERTIES  OF  ATMOSPHERE  ALONG 

THE  FLIGHT  PATH  AND  AT  BARCROFT 

(2)  OTHER  METEOROLOGICAL  PARAMETERS  AT  BARCROFT 


BACKGROUND 


I  ()N(;  JUMP  i 

•  METEOROLOGICAL  DATA 

•  NEPHEI.ON^  ETER  -  VISmiLITY 

•  IMMEDIATE  BARCROFT  AREA 
LONG  JUMP  II 

•  METEOROLOGICAL  DATA 

•  NEPHELOMETER  VISIBILITY 

•  AEROSOL  SIZE  DISTRIBUTION 

•  IMMEDIATE  BARCROFT  AREA 

•  MAY  NOT  REPRESENT  LINE  OF  SIGHT 

•  OROGRAPHIC  EFFECTS 

•  DIFFERENCES  IN  AIR  PARCELS 

LONG  JUMP  III 

•  AN  INSTRUMENTATED  AIRCRAFT  WAS  USED  TO 

SAMPLE  ALONG  OR  NEAR  THE  FLIGHT  PATH 


MEASUREMENTS 


•  AT  BARCROFT  LABORATORY,  THE  FOLLOWING 

PARAMETERS  WERE  MEASURED 

•  AEROSOL  PARTICLE  SIZE  DISTRIBUTION 

•  PARTICLE  SCATTERING  COEFFICIENT  (NEPHELOMETER) 

•  RELATIVE  HUMIDITY  (WET  AND  DRY  BULB  TEMPERATURES) 

•  AIR  TEMPERATURE  VERSUS  SOIL  SURFACE  TEMPERATURE 

•  INSOLATION 

•  WIND  SPEED  AND  DIRECTION 

•  OVER  BARCROFT  AND  ALONG  THE  FLIGHT  PATH,  THE 

FOLLOWING  PARAMETERS  WERE  MEASURED; 

•  AEROSOL  PARI  ICLE  SIZE  DISTRIBUTION 

•  AEROSOL  CHEMICAL  COMPOSITION 

•  AIR  TE M PE R ATUR E  AND  DEW  POINT 

•  CARBON  DIOXIDE  CONTENT 


CAL* 


CALCULATION  OF  Si 
COEFFICIENTS  FOR  IF 
WAVELENGTHS 

•  MIE  CODE  FROM ’F  ' 

.  INPUT 

•  RELATIVE'  :v 
.  PARTiaESEa 

•  REFRACTIVE  ii'TDICcS 


Exnnction.  1/km 


ion,  1/km 


Possible  causes  of  llucliiati>ms  at  e 

1.  Noise.  APS33  inlegraliuu  (.rr  tj,.:  <  pi.iiVi  2  inin.  :in<! 

that  of  the  EAA  was  i  iiiiii.,  l>i;t  j><  i  ■•■liicit-.  <•=  il;  (!:ita  is  aldii* 

6  min. 

2.  Terrain  contour  ellect.s.  Aerosols  lo.iy  !>■•  v  lir.i-iiii  ain!  in  tin 

vicinity  of  7000  ft  peaks  alonj;  nigi.t  and  (koktec! 
elsewhere. 

3.  Convection  cells.  Cells  can  i)e  c\[ocie(-  ^  5  ti,  is  ,)ni  v.  idc 

over  the  desert,  during  sjniu’ie”';.!!.  o  •,  j 2,,;, 
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Accumulalion  Mode  Aerosol  Conceiilralion  from  F.AA  I-'liglit  1 5  ai  1118  to  1 238  F’DT  on  25  August 


Outbound  Distance  (nm)  from  White  Mountain  Return 


Fluctuation  Widths 


Week  1.  The  average  fluctuation  width  was  26.6  nm,  whicli  was 
determined  from  the  extinction  minima.  The  lluctuation 
widths  were  determined  from  tlie  APS33  data. 

Week  2,  The  average  fluctuation  width  was  19.4  nm  determined 
from  the  particle  concentration  minima  of  the  higher  restdiition 
FA  A  data. 


Altitude,  ft 


A  METHOD  OF  ESTIMATING  SURFACE  METEOROLOGICAL 
RANGES  FROM  SATELLITE  AEROSOL  OPTICAL  DEPTHS 

D.R.  Longtin’,  E.P.  Shcttlc^*,  and  J.R.  Hummel' 

('SPARTA,  Inc.,  24  Hartwell  Avenue,  Lexington,  MA  02173,^Gcophysics 
Laboratory/OPA.  Bedford,  MA  01731) 

Measurements  of  aerosol  optical  depth  from  NOAA  Advanced  Very  High  Resolution  Radiometers 
(AVHRR)  aboard  polar-orbiting  satellites  have  been  used  to  estimate  die  surface  meteorological  range 
over  oceans.  To  do  this,  a  lookup  table  of  optical  deplli  versus  meteorological  range  has  been  developed 
which  is  based  on  the  aerosol  extinction  profiles  in  LOWTRAN  7.  To  better  simulate  conditions  over 
occaas,  the  boundary  layer  height  has  been  lowered  to  0.5  km  and  an  extinction  profile  corresponding 
to  a  150  km  meteorological  range  has  been  added. 

The  lookup  table  has  been  applied  to  AVHRR  aerosol  optical  depths  near  .selected  island  locations 
where,  in  turn,  the  inferred  meteorological  ranges  were  validated  agaiast  surface  observations.  Results 
show  reasonable  success  in  locating  regions  where  meteorological  ranges  exceed  30  km,  but  more  precise 
classification  may  not  be  possible  due  to  uncertainties  in  both  die  measurements  and  observations. 
Unfortunately,  the  algorithm  docs  not  fare  as  well  when  the  observed  meteorological  ranges  arc  less 
than  30  km. 

‘Now  at  NRL. 


A  METHOD  OF  ESTIMATING  SURFACE  METEOROLOGICAL  RANGES 
FROM  SATELLITE  AEROSOL  OPTICAL  DEPTHS 

Presented  at 

The  Annual  Review  Conference  on  Atmospheric  Transmission  Models 

Juno  6.  1990 

David  R.  Longtin^  Eric  P.  Shettie-  and  John  R.  Hummel' 

‘SPARTA.  Inc. 

24  Hartwell  Avenue 
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PURPOSE  OF  RESEARCH 


•  Determine  Whether  Measurements  of  Aero.sol  Optical  Depth  From  Polar- 
Orbiting  Satellites  Can  Be  Used  to  S.c'iir.iX  Sodace  Meteorological  Ranges 
Over  Oceans 


These  Are  Preliminary  Results 


METHODOLOGY 


•  Obtain  Experimental  Global  Sets  of  Aerosol  Optical  Depth  Over  Oceans  From 
NOAA/NESDIS 

•  Develop  a  Lookup  Table  of  Surface  Meteorological  Range  Versus  Aerosol 
Optical  Depth 

•  Extract  Aerosol  Optical  Depths  Near  Islands  and  Convert  Them  to  Surface 
Meteorological  Ranges 

•  Validate  Inferred  Surface  Meteorological  Ranges  with  Observed  Values 


DESCRIPTION  OF 
AEROSOL  OPTICAL  DEPTHS 


•  Based  on  Radiance  Measurements  From  the  Advanced  Very  High  Resolution 
Radiometer  (AVHRR)  Aboard  NOAA  Sun-Synchronous  Polar-Orbiting 
Satellites 

•  Retrieval  Algorithm  Relates  Aerosol  Optical  Depth  to  the  Upwelling  Radiation 
in  Channel  1  (0.58/im-0.68/im) 

•  Measurements  Are  Taken  When  NOAA  Satellites  Cross  the  Equator  (About 
13:30  Local  Time) 


Reported  Aerosol  Optical  Depths  Are  Scaled  to  0.50  /an 


Note:  Profiles  Based  on  Modified  Shetlle  and  Fenn  (1976) 

Note:  Profiles  Are  For  Surface  Meteorological  Ranges  of  150,  50,  23,  10  and  5  km 
Note:  Surface  Meteorological  Range  Equals  3.912/Totat  Extinction  at  0.55 /^m 


Note:  Dashed  Line'^  R-cablian  .1  H  ;  -  '  m,.  ■  :  !-;  '  AVHRR 

Measurement'}  (  •  rj  05,; 


VALIDATION  SCHEME 


SfAftTAWC 


•  Compute  Daily  Averages  of  Aerosol  Optical  Depth  Near  Island  Locations 

-  Use  All  Available  Data  Within  5  Degrees  of  Each  Island 

•  Omit  Those  Daily  Averages  Having  Excessive  Variability 

•  Convert  Average  Optical  Depths  to  Inferred  Meteorological  Ranges  Using 
Lookup  Table 

•  Obtain  Observed  Meteorological  Ranges  For  Days  When  Inferred  Values  Are 
Available 

-  Based  on  Observations  of  Surface  Visibility  Provided  by  USAFETAC 

•  Determine  the  Number  of  Times  When  the  Inferred  and  Observed 
Meteorological  Ranges  Are  in  Agreement 


SMRTAMC 


ISLAND  LOCATIONS  USED  FOR 
VALIDATION  STUDIES 


Ascension  Island 


•  a  3  c 


^  HISTOGRAMS  OF  INFERRED  AND  OBSERVED 
METEOROLOGICAL  RANGE 

SPARTA  INC. 
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FACTORS  LIMITING  VALIDATION  FOR 
METEOROLOGICAL  RANGES  >  30  KM 


•  Island  Locations  Often  I  ack  Reference  Points  for  Very  Clear  Conditions 


•  Maximum  Reported  Metecr.,!ioQ'oa;  Range  !s  Station  Dependent 

-  Marshall  Islands:  31  km 

-  Honolulu,  HI;  52  km 

-  Cantania.  Sicily:  39  ko! 

-  Jamaica:  30  km 

-  San  Juan,  PR:  46  km 

-  Tahiti;  52  km 

-  Ascension  island:  97  km 

•  For  Very  Clear  Cond'tions.  Av-raoe  Opt:"  k  Dopttis  Are  Similar  In  Magnitude 
to  the  Uncertainty  in  iris'  AVriRrl  As  io.io’  Optic.ii  Ds^pthc  (  •  C.05) 


CONCLUSIONS 


► 

SIVkflTAMC 


Reasonable  Success  in  Locating  Regions  Where  Meteorological  Ranges 
Exceed  30  km 


-  More  Precise  Classification  May  Not  Be  Possible  Due  To  Uncertainties  in 
Both  the  Measurements  and  Observations 


Limited  Success  in  Locating  Regions  of  Low  Meteorological  Range 

-  High  Meteorological  Ranges  Are  Often  Inferred  When  the  Obsen/ed  Value 
is  Below  20  km,  Possibly  Due  To  Inhomogeneties  in  Aerosol  Loading 

-  Improved  Results  Could  Be  Expected  If  It  Were  Possible  to  Define  the  Top 
of  the  Boundary  Layer  More  Accurately 


A  COMPARISON  OF  THE  UVTRAN  AND  LOWTRAN  7  AEROSOL  MODELS 

S.  O’Brien 

Las  Cmces  Scientific  Consulting,  Las  Cnices,  NM  88001 

A  comparison  between  aerosol  model  results  used  by  the  AFGL  LOWTRAN  7  and  ASL  UVTRAN 
models  is  examined  for  a  variety  of  surface  visibilities.  Recent  semicmpirical  results  from  particle  sizing 
problems  operated  near  the  surface  at  the  White  Sands  Missile  Range  are  also  compared  with  the  model 
results. 


A  Comparison  of  the  UVTRAN  and  LOWTRAN  7 

Aerosol  Models 


Sean  G.  O’Brien 

Las  Cruces  Scientific  Consulting 
6  June  1990 


Summary  of  Presentation 


•  Description  of  UVTRAN  aerosol  model 

•  Description  of  LOWTRAN  7  aerosol  models 

•  Review  of  past  UVTRAN  -  LOWTRAN  aerosol  comparisons 

•  Comparison  of  UVTRAN  and  LOWTRAN  7  for  sea  level  aerosols 

•  Comparison  of  UVTRAN  and  LOWTRAN  7  for  elevated  desert 
aerosols 

•  Comparison  of  model  results  with  measurements  made  in  an 
elevated  desert  region 


•  Conclusions 


UVTRAN  Aerosol  Model 


•  References:  E.M.  Patterson,  'Ultra-Violet  Atmospheric 

Propagation  Model,'  Final  Report  on  Project 
DAAL03-86-D-001  Delivery  Order  0578  (1988) 

E.M.  Patterson  and  J.B.  Gillespie,  Appl.  Opt.  28, 
p.  245  (1989) 

•  Aerosol  model  spans  185  to  700  nm  wavelength  region 

•  Aerosol  extinction  is  dominant  In  the  300  to  700  nm  region 


where  q  -  0.585  V 


LOWTRAN  7  Aerosol  Models 


•  References:  F.X.  Kneizys,  et  al..  LOWTRAN  6  technical 

report.  AFGL-TR-83-0187  (1983) 

F.X.  Kneizys.  et  al..  ’User's  Guide  to 
LOWTRAN  7.’  AFGL-TR-88-0177  (1988) 

D.R.  Longtin.  E.P.  Shettle.  J.R.  Hummel,  and 
J.D.  Pryce,  *A  Wind  Dependent  Desert  Aerosol 
Model:  Radiative  Properties,'  AFGL-TR-88-01 12 


•  Models  considered  here: 

Sea  Level  Elevated  Desert 


Urban,  70%  RH  Rural,  0%  RH 

Maritime,  70%  RH  Desert,  WS-0  m/s 

Rural.  70%  RH 


A  Review  of  UVTRAN  -  LOWTRAN  6 
Aerosol  Comparisons 


•From  Patterson  and  Gillespie,  1989 


•  Conditions: 


•  Results: 


AFGL  Rural  aerosol  model,  70%  RH 

P  -  1013  mb.  T  ■  2880K 

V  -  23  km 
560 


Wavelength  k  (km  ^ ) 

(nm) 


UVTRAN 

LOWTRAN  6 

300 

0.434 

0.276 

550 

0.158 

0.158 

700 

0.106 

0.1  18 

UVTRAN  -  LOWTRAN  6  Comparison 

(continued) 


Reference:  W.A.  Baum  and  L.  Dunkelman,  'Horizontal 

Attenuation  of  Ultraviolet  Light  by  the  Lower 
Atmosphere,'  JOSA  45,  166  (1955) 

Results  shown  here  are  from  Patterson  and  Gillespie,  1969 
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1.21 
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Comparison  of  UVTRAN  with  LOWTRAN  7 
Sea  Level  Aerosols 

Conditions:  P  ■  1013  mb.  T  -  288®K.  RH  -  70% 


Results  for  300  nm: 
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Comparison  of  UVTRAN  with  LOWTRAN  7 
Sea  Level  Aerosols 


Results  for  700  nm: 


550 


k  (km"') 
ext 


(km) 

Urban 
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Rural 

UVTRAN 
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0.144 

0.1 18 
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0.0393 

Comparison  of  UVTRAN  with  LOWTRAN  7 
Elevated  Desert  Aerosols 


Conditions: 

P  ■  880  mb,  T 

•  303  °K,  RH 

o 

■ 

Results  for 

300  nm: 

k  (km  S 

V 

550 

ext 

(km) 

Rural 

Desert 

UVTRAN 

60 

0.120 

0.122 

0.254 

100 

0.0517 

0.0527 

0.154 

150 

0.0290 

0.0295 

0.109 

Comparison  of  UVTRAN  with  LOWTRAN  7 
Elevated  Desert  Aerosols 


Results  for  700  nm: 
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_(Lrni_ 

Rural 

Desert 

UVTRAN 
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0.0514 

0.0493 

0.0409 

100 

0.0221 

0.02  1 2 

0.0154 

150 

0.0124 

0.0119 

0.00785 

Comparison  of  Model  Results  with 
Elevated  Desert  Measurements 

References;  J.  Boatman,  D.  Wellman,  and  B.  Bodhaine  of 
NOAA/WPL,  D.  Garvey  of  ASL,  private 
communications,  1989 

S.  O'Brien  and  D.  Longtin,  OMI-371,  1989 

Measurements  examined  here  are  of  two  kinds: 

-  Aircraft  three  color  nephelometer  that  provided 
aerosol  scattering  at  450,  550,  and  700  nm 

-  Surface  Knollenberg  particle  sizing  probes  (LAS-X 
type)  that  provided  size  distribution  inputs  for 
Mie  computations 


Dn/Dr  (cm'^  /irr, 


1A;:-X  Aircrufl  Duta, 


Example  of  least-squares  bimodal  lognormal 
curve  fit  with  automatic  edit 


Seasonal  Size  Distribution  I’araoctcrs  Used  in 
boundary  Layer  Kie  Calculations 
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Conclusions 


For  sea  level  aerosols,  visibilities  greater  than  10  km, 
and  short  wavelengths  (<  550  nm),  UVTRAN  predicts  a  higher 
aerosol  extinction  than  LOWTRAN  7;  limited  experimental 
data  tend  to  support  UVTRAN  in  such  conditions 

For  sea  level  aerosols,  the  two  models  agree  reasonably 
well  when  either  the  visibility  is  low  or  the  wavelength 
is  long 

For  elevated  desert  aerosols,  the  two  models  are  in 
agreement  only  when  both  the  visibility  is  relatively  low 
(<  23  km)  and  the  wavelength  is  above  550  nm 

Under  the  high  visibility  conditions  of  elevated  deserts, 
the  UVTRAN  model  predicts  significantly  stronger  spectral 
dependence  of  aerosol  extinction  than  does  LOWTRAN  7; 
limited  experimental  data  indicate  that  the  LOWTRAN  7 
desert  aerosol  model  is  more  accurate  in  such  conditions 


HIGH  RESOLUTION  SOLAR  SPECTRUM  BETWEEN  2000  AND  3100  ANGSTROMS 

L.A.  Hall  and  G.P.  Anderson 

Geophysics  LaboratoryA-IM/OPE,  Hanscom  AFB,  MA  01731 

Solar  spectra  in  the  wavelength  range  2000-3100  Angstroms,  measured  at  40  km  in  the  stratosphere, 
have  been  extrapolated  to  zero  optical  depth  to  provide  a  reference  spectrum  in  0.1  Angstrom  resolution. 
The  wavelength  scale  has  been  carefully  compared  to  wavelength  standards  and  a  cyclic  instrumental 
effect  removed,  resulting  in  a  wavelength  accuracy  to  within  .04  Angstroms.  The  absolute  intensities 
have  been  normalized  to  a  previously  released  spectrum  in  1.0  Angstrom  resolution. 


HIGH  RESOLUTION  SOUR  SPECTRUM  BETWEEN  2000 
AND  3100  ANGSTROMS 


L.  A.  HALL  AND  G.  P.  ANDERSON 
GEOPHYSICS  UBORATORY 


Balloon  Measureirients  of  Stratospheric  Ultraviolet 


0  Five  Flights:  21  April,  1977 

19  April,  1976 
27  April ,  1960 
22  April.  1981 

20  April,  1983 


0  Results:  Absolut'’  irradiance  in  2000-3100  A  range 

at  balloon  altitudes 

Variation  of  UV  irradiance  in  11-yr  cycle 

In  situ  measurement  of  Heraberg  continuum 
cross  section 

Measurement  of  02  transmission  in  the 
stratosphere 

Photodissociati on  rate  coefficients  of  02 
in  the  stratosphere 

Extrapol.‘tion  of  absolute  irradiance  to 
zero  optical  depth 
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PREDISSOCIATION  LINE  WIDTHS  OF  THE  SCHUMANN-RUNGE 
ABSORPTION  BANDS  OF  SOME  ISOTOPES  OF  OXYGEN 

W.H.  Parkinson',  A.S.-C.  Cheung^,  S.L.  Chiu^,  J.R.  Esmond',  D.E.  Freeman  and  K.  Yoshino' 
'Harvard-Smithsonian  Center  for  Astrophysics,  Cambridge,  MA  02138, 

^Chemistry  Department,  Univeisity  of  Hong  Kong,  Hong  Kong 

Predissociation  line  widths  of  the  (3,0)  -  (1 1,0)  Schumann-Runge  bands  of  '*©2  and  '^0'*0  in  the 
wavelength  region  189-196  nm  have  been  obtained  from  the  published  measurements  of  the  absolute 
absorption  cross  sections  of  Yoshino  et  al.  [Planet.  Space  Sci.,  36,  1201  (1988);  37,  419  (1989)]  and 
spectroscopic  constants  of  these  molecules  of  Cheung  [J.  Mol.  Spectrosc.,  131,  96  (1988);  134,  362 
(1989)].  The  line  widths  arc  dctcmiined  as  parameters  in  the  non-linear  least  squares  fitting  of  calculated 
to  measured  cross  sections.  Predissociation  maxima  arc  found  at  upper  vibrational  levels  with  v’=  4,  7  and 
10  for  '*02  and  '*0'*0.  Our  predissociation  line  widths  arc  mostly  greater  than  previous  experimental 
values  for  both  isotopic  molecules.  This  work  is  supported  by  NASA  grant  NAG  5-484  to  Smithsonian 
Astrophysical  Observatory. 


Prcdissociatiun  Line  Widths  of  tlie  Schuinann-Runge  Absorption  Dands 
of  Some  fsotopcs  of  Oxygen 

W  H.  Parkinson',  A  S  C.  Cheung*,  S  L.  Chiu*.  JK  tZsniond',  D.K.  Freeman'  and 
K.  Yoshino' 

'Harv&rd-Smithsonixn  Center  for  Asirophyaica.  Cambridge,  MA  02138 
’Chemistry  Department,  University  of  Hong  Kong,  Hong  Kong 


Prediaaociation  line  widths  of  the  (3,0)-(n,0)  Schumann-Runge  bands  of  ‘*0)  and 
‘*0**0  in  the  wavelength  region  180-196  nm  have  been  obtained  from  the  published 
meaaurcmenta  of  the  absolute  absorption  cross  sections  of  Yoshino  ^  [PlaiMt  Sf>ace  Sci 
1201  (1988);  419  (1989)}  arvd  spectroscopic  constants  of  these  molecules  of 

Cheung  |J.  Mot.  S^>cctrasc  131 .  96  (1988),  134,  362  (1989)]  The  line  widths  are 

determined  as  parameters  in  the  non  linear  least  squares  fitting  of  calculated  to  measured 
cross  sections  Prcdissociatiun  maxima  are  found  at  upper  vibrational  levels  with  v'b  4,  7 
and  10  for  '*0]  and  '*0'*0.  Our  predissocialion  line  widths  are  mostly  greater  than 
previous  experimental  values  for  both  isotopic  molecules.  This  work  is  supported  by 
NASA  grant  NAG  5-484  to  Smithsonian  Asirophysical  Observatory 


Predissociation  Line  Widths  of  the  Schumann-Runge  Absorption  Bands 
of  Some  Isotopes  of  Oxygen 


W.H.  Parkinson,  A.S-C.  Cheung,*.  S.S-L.  Chiu,*,  J.R.  Esmond, 
D.E.  Freeman  and  K.  Yoshino 


Harvard-Smi thsoni an  C<'nter  for  Astrophysics 
60  Garden  Street,  Caml)ridge.  MA  02138,  USA. 

‘Chemistry  Department,  University  of  Hong  Kong,  Hong  Kong 


Supported  by  the  NASA  Upper  Atmo-spherlr  Re-searcli  Prograa 


I’hf  Scliiiiuaim  Rung!'  ahsorplicju  hands  of  **’02  and  their  extensive 
predi,ss<Kiation  are  of  considerable  atmospheric  significance  in  connection  with  the 
triuismission  of  solar  radiation  and  the  production  of  0(^P)  atoms.  The  cross  sections 
of  the  S-R  bands  consists  of  hundreds  of  broadened  and  often  overlapping  rovibronic 
lines,  and  the  predissociation  of  “'O2  is  the  main  source  of  odd  oxygen  in  the 
atmosphere  above  GO  km.  The  cross  sections  are  also  temperature  dependent  and  can 
only  be  obtained  from  {xirameters  including  predissociation  line  widths  at  any 
temperature.  The  predissociation  line  widths  and,  in  particular,  their  variation  with 
the  up{>er  level  vibrational  quantum  number  v',  are  also  of  importance  in  the 
elucidation  of  the  curve-crossing  mechanism  resfionsible  for  the  predissociation.  The 
analogous  predissociation  line  widths  of  the  Scliumann-Runge  bands  of  *®02  and  *®0*®0 
would  provide  valuable  supplementary  information  on  the  meclianism  of 

predissociation  because  the  predissociating  vibronic  energy  levels,  but  not  the  associated 
Ixxind  and  repulsive  potential  energy  curves,  are  isotope-dependent.  In  addition, 
constitutes  about  0.4*’x)  of  atmospheric  O2  and  is  a  source,  via  predissociation,  of 
**()  in  the  .stratosphere. 


THE  PREDISSOCIATION  OF  B  STATE  OF  Oj. 

The  predi-s-sociation  i.^  caused  by  four  repulsive  states,  namely,  *11^,  ^Ily,  ^II^  and 
which  cross  the  upper  I)  state  of  the  Sdiumann-Runge  system  in  its  bound 
region  The  theoretical  investigations  of  Julienne  and  Krauss  (1975)  and  Julienne 
(I'JTGJ  have  shown  th.it  both  the  level  shifts  and  line  widths  can  be  attributed  to 
(lert urb.it ions  dominated  by  the  '^11^  state,  with  other  states  playing  minor  roles  in  the 
predi.s'ix  iation.  Nevertheless,  accurate  predis-soi  iation  line  width  metisuremeiiLs  for 
the  v.irioti,^  vibration.il  levels  are  required  to  provide  a  meimingful  coiniiiirison  of 
ililferent  .-'et.s  of  line  widths  calculated  from  theoreticid  piirameters. 
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JULIENNE  AND  KRAUSS 
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PARAMETERIZATION  OF  THE  SCHUMANN-RUNGE  BANDS  (Absorption  Cross  Section) 


1.  Boltzfflan  Population 

2.  Honl  London  Factor 

3.  Band  Oscillator  Strengths 

‘‘Oj;  Planet.  Space  Sci  35,  1067-1075  (1987). 

‘*0j;  Planet.  Space  Sci.  36,  1201-1210  (1988). 

Planet.  Space  Sci  37,  .119-426  (1989). 

1.  Line  Center  Positions 
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5.  Predissociation  Line  Widths 

'"Oj;  ,1  Chem  Phys  TJ ,  812-819  (199"; 

'“Oj;  Present  work 
'^O^O,  Present  work 
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TABLE  1.  Predissociation  line  widths  (cm'*,  FWHM)  of  the  Scliumaiin-Runge 
bands  of  *®02  and  **^0**0 


“Oj 

'»Oj 

lOQiso 

1 

0.93±0.18 

2 

0.72±0  04 

3 

1,78±0.11 

1.53±0.15 

1.62±0.20 

4 

3.87±0  27 

3.23±0.11 

3.69±0.18 

5 

2.13±0.07 

3.12±0.08 

2.68±0.13 

6 

1.79±0.21 

1.05  ±0.03 

1.32±0.10 

7 

2,01±0.10 

2.77±0.11 

2.47±0.20 

8 

1.92±0  12 

1  32±0.18 

1.62±0.20 

9 

1,01±9.16 

1.1 8  ±0.07 

1.00±0.07 

10 

1.09±0.08 

1.71  ±0.07 

1.67±0.10 

11 

1.48±0  18 

0.84  ±0.05 

1.34±0.18 

12 

0.88  ±0.1 4 

Variation  of  pr  e  J  i  s  soc  »  a  t  1  on  lln*  »idlh  wit»  v  i  tr  a  ^  i-'na  I 
nijfflber  In  the  8  ’Ej  »tal«  of  ’*0,  FjIltfJ  c’.tcl^s  prestn  a.>*:w»- 
values,  op*»n  circles:  Lewis  sol  l>1  Ilr.r  Julienne 


AURIC:  ATMOSPHERIC  ULTRAVIOLET  RADIANCE  INTEGRATED  CODE 

R.L.  Huguenin,  M.A.  LeCompie 

Aerodyne  Research,  Inc,,  45  Manning  Road,  Billerica,  MA  01821 

R.E.  Huffman 

Geophysics  Laboratory/LIM,  Hanscom  AFB,  MA  01731 

AURIC  is  a  planned  computer  simulation  model  to  synthesize  background  radiance  and  luttcr 
as  viewed  by  ultraviolet  sensors.  Radiometrically  “accurate”  radiance  and  clutter  predictions  will  be 
emphasized  through  inclusion  of  all  relevant  phenomenology.  Radiances  will  be  fully  traceable  to 
underlying  phenomenology  using  parametric  first  principles  radiance  models.  Included  will  be  the 
capability  to  compare  AURIC  predictions  with  measured  data.  Phenomenological  models  and  simulation 
algorithms  will  be  modularized  to  accommodate  upgrades  and  substitutions.  A  relational  data  base 
management  system  will  serve  as  the  shell  for  integration  and  control  of  the  modules,  as  well  as  provide 
the  foundation  for  a  high  level  user  interface  that  provide*^  case  of  use  for  scientists  and  system  engineers. 


AERODYNE  RESEARCH,  Inc. 

SYSTEMS 

AURIC:  ATMOSPHERIC  ULTRAVIOLET  RADIANCE  INTEGRATED 

CODE 


■  R.L.  HUGUENIN,  ‘M.A.  LeCOMPTE,  AND  “R.E.  HUFFMAN 


•AERODYNE  RESEARCH,  INC. 
45  MANNING  ROAD 
BILLERICA,  MA  01821 


•  •  IONOSPHERIC  MODELING  AND  REMOTE  SENSING  BRANCH 
IONOSPHERIC  RESEARCH  DIVISION 
AIR  FORCE  GEOPHYSICS  LABORATORY 
HANSCOM  AFB,  MA  01731 
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SYSTEMS 


PROBLEM  BEING  ADDRESSED 

•  MEASUREMENTS  LIMITED 

-  TEMPORAL  COVERAGE 

-  SPATIAL  COVERAGE 

•  SENSOR  CHARACTERISTICS 

•  MEASUREMENTS  OF  ALL  BACKGROUND  CONDITIONS  UNREALISTIC 

•  RADIANCE  MODELING  CAPABILITY  NEEDED 

INTERE'OLATE  AND  EXTRAPOLATE  ML  ASUREO  RADIANCES  TO 
OTHER  CONDITIONS 

-  SYNTHESIZE  RADIANCES  FOR  SCENARIOS  SENSORS  FOR  WHICH 
NO  DATA  EXIST 

ANALYZE  AND  TEST  CURRENT  UNDERSTANDING  OF  SOURCES 
OF  ULTRAVIOLET  BACKGROUND  RADIANCE  AND  CLUTTER 


PROVIDE  INTEGRATED  SVNTHETIC  SCENES,  WHERE  RADIANCE  INTERACTIONS 
ARE  TOTAL!  Y  CONSISTENT  AND  RADIOMETRICALLY  ACCURATE 


liYSTEl/iS 


AURIC 


UV  FIELD 
MEASUREMENTS 


UV  SENSOR 
SYSTEMS  ANALYSIS 


UV  LABORATORY 
AND  SIMULATION 
MEASUREMENTS 
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SYSTEMS 


AURIC  PROGRAM  OBJECTIVE 


DEVELOP  AND  PROVIDE  TO  AIR  FORCE  GEOPHYSICS  LABORATORY 
A  COMPUTER  SIMULATION  MODEL  TO  SYNTHESIZE  BACKGROUND 
RADIANCE  AND  CLUTTER  AS  VIEWED  BY  USER  SPECIFIED  ULTRA¬ 
VIOLET  SENSORS 


SYSTEMS 


AURIC  REQUIREMENTS  (CONT"D) 


•  HIERARCHICAL  DEVELOPMENT  STRATEGY 

-  NEED  TO  ACCOMMODATE  FUNDING  LEVEUSCHEDULE 
UNCERTAINTIES 

-  COMPLETE  “ACCURATE"  MODEL  FOR  RESTRICTED  SETS  OF 
MEASUREMENT  SCENARIOS 

-  ADD  MEASUREMENT  SCENARIO  VERSATILITY  ACCORDING  TO 
FUNDING  AVAILABILITY 

•  ACCESSIBLE  DESK-TOP  CAPABILITY 

-  ANALOGOUS  TO  PC  LOWTRAN.  PC  SPIRITS 

-  EASE  OF  USE  FOR  SCIENTISTS,  SYSTEMS  ENGINEERS 

•  FULLY  SUPPORTED,  DOCUMENTED 

•  RESPONSIVE  BUT  CONTROLLED  UPDATE  RELEASESA/ERSIONS 

-  ENCOURAGE  USER  CONTRIBUTIONS 

-  CONTROLLED  RELEASES  TO  PERMIT  UNIFORM  COMPARISONS 

-  PERIODIC  USER  MEETINGS 


SYSTEMS 


^1- 


AURIC  REQUIREMENTS 

•  RAOIOMETRICALLY  "ACCURATE”  RADIANCE  AND  CLUTTER  PREDICTIONS 

ALL  RELEVANT  PHENOMENOLOGY  NEEDED  FOR  USER-SPECIFIED 
SCENARIOS 

ENTRANCE-APERTURE  SPECTRAL  RADIANCES  NEED  TO  EXCEED 
SENSOR  SPECIFICATIONS 

CAPABILITY  NEEDED  TO  SPECIFY  SENSOR  MODEL  AND  APPLY  TO 
RADIANCE  PREDICTION 

INCLUSION  OF  LOCAL  ENVIRONMENT  CONTAMINANTS  DESIRED 

•  RADIANCE  FULLY  T  HACEABLE  TO  UNDERLYING  PHENOMENOLOGY 

REQUIRES  FULLY  PARAMETRIC  FIRST-PRINCIPLES  RADIANCE 
MODELS 

AB  INITIO  CALCULATIONS 
•  PARAMETRIC  PERTURBATION 

•  CAPABILITY  TO  COMPARE  AURIC  PREDICTIONS  TO  MEASURED  DATA 


REQUIREG  CAPABILH  /  TO  DISPLAY  BOTH  DATA  SETS 
HEQUIREG  CAPABILITY  TO  COMBINE  BOTH  DATA  SETS 
VALiDATION  CAl  IBFIATION  MODEL  REFINEMENT  CAPABILITY 
DESIRED 
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PHENOMENOLOGICAL 

MODELS 


SIMULATION 

ALGORITHMS 


PHYSICAL 

PROPERTIES 

MODULE 


PARAMETER 

SPECinCATTON 

MODULE 


RADIANCE 

MODULE 


/  RADIANCE  \ 
\  TABLES  I 


I  PARAMETERS  ^ 
\  DATA  A 


3D  SPATIAL 
SYNTHESIS 
MODULE 


f  SPATIAL  \ 
%  DESCRIPTION  i 


RADIANCE 

SIMULATION 

MODULE 
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SPECTRAL  RADIANCE 
OUTPUT  DATA 


SYSTEMS 


AURIC  SYSTEM  DESIGN/DEVELOPMENT  APPROACH 

SEPARABLE  MODULAR  DESIGN 

-  PHENOMENOLOGY  UPDATBSUBSTITUTION  ACCOMMODATION 

-  HIERARCHICAL  DEVELOPMENT  ACCOMMODATION 

-  EXPANDABILITY 

-  SHARED  DEVELOPMENT  EFFICIENCY 
SIMULATION  MODELING 

-  COMPUTATIONAL  EFFICIENCY  REQUIREMENT 

-  RADIOMETRIC  ACCURACY  NOT  MEASURABLY  COMPROMISED 

-  PHENOMENOLOGICAL  TRACEABILITY  NOT  MEASURABLY 

COMPROMISED 

RULE-BASED  SCENE  SPECIFICATION  AND  MODELING 

-  ACCOMMODATES  COMPLEX  PHENOMENOLOGY  CONSISTENCY 

REQUIREMENTS 


PERMITS  HIGH  LEVEL  SCENE  SPECIFICATION  WITH  OPTIONAL 
I  OWER  LEVEL  ACCESS 


SYSTEMS 

PHYSICAL  PROPERTIES  MODULE  TABLE  ENTRIES 

•  RELEVANT  CROSS-SECTIONS  WITH  SPECTRAL  DEPENDENCE 
FOR  EACH  ATMOSPHERIC  SPECIES 

PHOTOABSORPTION 

PHOTODISSOCIATION 

PHOTOIONIZATION 

ELECTRON  IMPACT 

ION  IMPACT 

NEUTRAL  IMPACT 

METEORITIC  IMPACT 

COSMIC  RAY  (SOLAR,  GALACTIC)  IMPACT 

•  BRANCh.NG  RATIOS  FOR  RELEVANT  ATMOSPHERIC  SPECIES 

•  OTHER 


SYSTEMS 


AURIC  CODE  MODULES 


•  PHYSICAL  PROPERTIES  MODULE 

-  SETS  OF  LOOK-UP  TABLES  FOR  RELEVANT 
PHYSICAL  PROPERTIES 

-  LOOK-UP  INTERPOLATORS 

-  SELECTION  RULES  RESIDE  IN  PARAMETER 
SPECIFICATION  CODE 


RADIANCE  MODULE  (CONTINUED) 


RADIANCE  LOOK-UP  TABLES 

•  COMPUTATIONALLY  MOST  EFFICIENT  RADIANCE  SOURCE 

•  PREDETERMINED  SPECTRAL  RADIANCES  VS.KEY  PHYSICAL 
PARAMETERS 

-  COMPUTATIONALLY  INTENSIVE  CALCULATIONS  SHIFTED 
TO  OFF  -  LINE  SYSTEM 

-  MEASURED/CALIBRATED  SPECTRAL  RADIANCE  DATA  FOR 
KNOWN  CONDITIONS 

-  CALIBRATED  TO  COMMON  RADIANCE  UNITS 

-  REVISIONS/UPDATES  PROVIDED  AS  CONTROLLED  RELEASES 

•  PARAMETRIC  INTERPOLATION  EXPRESSIONS 

PARAMETRICALLY  SEPARABLE  APPROXIMATION  FORMULAS 

•  SELECTION  RULES/CONTROLLER  IN  PARAMETER 
SPECIFICATION  CODE 


I 


SYSTEMS 


AURIC  CODE  MODULES 

•  RADIANCE  MODULE 

-  RADIANCE  LOOK-UP  TABLES  GENERATED  BY  OFF  -  LINE  CODE 

-  RADIANCE  CALCULATION  CODES 

-  RADIANCE  PERTURBATION  CODES 


-  SELECTION  RULES/CONTROLLER  IN  PARAMETER 
SPECIFICATION  CODE 


SYSTEMS 


> 


3D  SPATIAL  SYNTHESIS  MODULE 

•  ANALYTICAL  SPECIFICATIONS  OF  MEAN  VALUE  FOR  EACH 
RELEVANT  PHYSICAL  PROPERTY 

LATITUDE,  LONGITUDE,  ALTITUDE  POINT  SPECIFICATION 

CONTINUOUS  RATHER  THAN  ARTIFICIALLY  LAYERED 

SPATIAL  RESOLUTION  FINER  THAN  SENSOR  IFOV 

•  SPATIAL  CLUTTER 

SPATIALLY  DISCREET  PHENOMENA  MODELED  WITH 
BOUNDED  TEXTURE  MODEL 

SPATIALLY  NON-DISCREET  PHENOMENA  PERTURBED 
FRACTALLY 

EMPIRICALLY/THEORETICALLY  DERIVED  FRACTAL 
TEXTURE  DIMENSION 

•  LINE-OF-SIGHT  INTERSECTION  OF  PHYSICAL  PROPERTIES 

VOXEUSCENE  ELEMENT  ANALYTICAL  INTERSECTION  MODEL 
RAY  TRACE  PHYSICAL  PROPERTIES  SAMPLING  ALGORITHM 
1/4  -  VOXEL  PHYSICAL  PROPERTIES  AVERAGING  ALGORITHM 


SYSTEMS 


RADIANCE  MODULE  (CONTINUED) 


•  RADIANCE  CALCULATION  CODES 

DIRECT  CALCULATION  FOR  COMPUTATIONALLY  EFFICIENT 
CODES 

•  RADIANCE  PERTURBATION  CODES 

-  PARAMETRIC  RADIANCE  EXPRESSIONS  FOR  EACH 
RELEVANT  PHENOMENON 

-  CLUTTER  RADIANCES  FROM  PHYSICAL  PROPERTY 
TEXTURES 


COMPUTATIONALLY  EFFICIENT 


SLANT  PATH  MANIPULATION 

IMAGE  PIXEL  INTERPOLATION  (FUTURE) 


PARAMETER  SPECIFICATION  MODULE 


•  USER  INTERFACE  TO  OVERALL  MODULE 

HIGH  LEVEL  SCIENTIFIC/ENGINEERING  SPECIFICATIONS 

OPTIONAL  LOWER  LEVEL  SPECIFICATIONS 
CONSTRAINED  USER  INPUTS  -  PHYSICAL  CONSISTENCY 

TRANSFORMATION  OF  USER  INPUTS  TO  CODE  INPUTS 

•  SYSTEM  CONTROLLER 

3D  SPATIAL  DISTRIBUTION  MODELS/RULES 

RADIANCE  MODELS/CALCULATION  RULES 
INTERNAL  CONSISTENCY  RELATIONSHIPS 
MODULE  DRIVER/SEQUENCE  RULES 

•  RELATIONAL  DATA  BASE  MANAGEMENT  SHELL 


SYSTEMS 


RADIANCE  SIMULATION  MODULE 


•  RADIANCE  ASSIGNMENT  TO  LINE-OF-SIGHT  PATH 

SEGMENTED  PATH  (VOXEL)  RADIANCE 
CALCULATIONS 

LINE-OF-SIGrIT  PATH  RADIANCE  PROPAGATION 

TWO  RADIANCE  FORMATS 

SPECTRAL  RADIANCE  AT  SENSOR  ENTRANCE 
APERTURE 

SPECTRAL  RADIANCE  CONVOLVED  WITH 
SENSOR  MULTISPECTRAL  TRANSFER  FUNCTION 


PARAMETER  SPECIFICATION  MODULE 


SYSTEMS 


•  3D  SPATIAL  DISTRIBUTION  MODEL  SOURCES 
MSIS  THERMOSPHERE  MODEL 
INTERNATIONAL  IONOSPHERE  MODEL 
US  STANDARD  ATMOSPHERE 
OTHER 


SYSTEMS 


PARAMETER  SPECIFICATION  MODULE  (CONT'D) 


•  3D  SPATIAL  DISTRIBUTION  MODELS 


NUMBER  DENSITY  DISTRIBUTION  MODELS  FOR  RELEVANT  ATMOSPHERIC 
SPECIES  (NEUTRAL,  CHARGED) 

NEUTRAL  ATMOSPHERIC  TEMPERATURE  DISTRIBUTION  MODEL 
ION  TEMPERATURE  DISTRIBUTION  MODELS 
PHOTOELECTRON  SPECTRUM  MODEL 
EXOSPHERIC  TEMPERATURE  MODEL 
SOLAR  CYCLE,  SOLAR  ACTIVITY  MODELS 
CORPUSCULAR  RADIATION  MODELS 


MAGNETIC  FIELD  MODEL 


SYSTEMS 


OUTPUT  AND  DISPLAY 

•  PARAMETER  SPECIFICATIONS 

INTERACTIVE  MENU,  PROMPTS,  FLAG, 

OUTPUT  SUMMARY 

•  SPECTRAL  RADIANCE  PLOT 

WITHOUT  SENSOR  FUNCTION  APPLIED 
WITH  SENSOR  FUNCTION  APPLIED 
SELECTABLE  RADIANCE  UNITS,  INCLUDING  SENSOR 
UNITS  (COUNTS,  DATA  NUMBER) 

SELECTABLE  WAVELENGTH/WAVENUMBER  UNITS 
NORMALIZED.  SCALED,  RATIO,  DIFFERENCE  OPTIONS 
SPECTRUM  SEQUENCE.  OVERLAY  OPTIONS 

•  CLUTTER  METRICS,  POWER  SPECTRAL  DENSITIES 

•  RASTER  IMAGE  DISPLAY  (FUTURE) 

•  IMAGE  MANIPULATION/ST ATISTICS/TRANSFORMS  (FUTURE) 


SYSTEMS 


PARAMETER  SPECIFICATION  MODULE 


RADIANCE  MODELS'CALCULATION  RULES 

•  SEASONAL  CYCLE  MODEL 

•  DIURNAL  CYCLE  MODEL 

•  PATH  ATTENUATION  MODELS  FOR  EACH  RELEVANT  SPECIES 

ABSORPTION 

SCATTERING 

•  PATH  EMISSION  MODELS  FOR  EACH  RELEVANT  SPECIES 

•  PATH  SCATTERED  RADIANCE  MODELS 


SINGLE/MULTIPLE  SCATTERING 


RESONANCE  SCATTERING 
•  OUTPUT  RADIANCE  INTERPOLATOR 


/  '1000 -4000  A  (100.000  -  25000  CM  )  INITIAL  RANGE 
5CM  '(0.05  -  0  80A)  INITIAL  SPECTRAL  RESOLUTION 


SYSTEMS 


AURIC  MODEL  DEVELOPMENT  SEQUENCE 


l-D(LINE-OF-SIGHT)  2-D  (IMAGERY) 


MID-LATITUDE  NADIR 

PHASE  I 

PHASE  11 

MID-LATITUDE  LIMB 

DISTURBED/AURORAL 

PHASE  III 

PHASE  V 

HIGH  SPECTRAL 

RESOLUTION 

PHASE  IV 

PHASE  VI 

SYSTEMS 


AURIC  CODE  DEVELOPMENT  ENVIRONMENT 

•  EVEREX  PC  (80386)  /MICROVAX  II 

•  OFF  -  THE  -  SHELF  RDBMS  (DEVELOPME(4T) 

ORACLE  RDBMS 

UNIX  WITH  X  -  WINDOWS/MOTIF 

FORTRAN  COMPILER 


•  PORTABILITY 


MEASUREMENT  NEEDS  FOR  AURIC:  ATMOSPHERIC 
ULTRAVIOLET  RADIANCE  INTEGRATED  CODE 

R.E  Huffman 

Geophysics  Laboratory.  Hanscom  AFB,  MA  01731 

The  development  of  the  \URIC  model  for  ultraviolet  radiance  and  transmission  will  require  the 
availability  of  a  wide  variety  of  both  laboratory  and  field  measurements.  These  measurements  or 
phenomenology  needs  are  available  for  many  of  the  areas  of  interest,  but  a  number  of  space  measurements 
have  not  been  conducted  over  a  sufficient  period  of  time  (solar  cycles)  and  with  well-calibrated  sensors. 
The  current  and  near-future  availability  of  airglow,  auroral,  and  scattering  space  measurements  of  use 
for  the  validation  and  improvement  of  the  AURIC  model  will  be  discussed.  Improved  measurements 
of  solar  radiation,  atmospheric  density  and  composition,  and  celestial  sources  will  all  contribute  to  an 
improved  model.  Laboratory  measurements  of  the  photon  cross  sections  for  the  Schumann- Runge  bands 
of  oxygen  that  have  recently  become  available  are  an  example  of  the  measurements  needed  for  the  model 
at  wavelengths  shorter  than  the  current  limit  of  LOWTRAN  7  of  200  nm. 


MEASUREMENT  NEEDS  FOR  AURIC: 
ATMOSPHERIC  ULTRAVIOLET  RADIANCE  INTEGRATED  CODE 


Annual  Review  Conference  on  Atmospheric  Models 
5-6  June  1990 
Geophysics  Laboratory 


Robert  E.  Huffman 
Geophysics  Laboratory 
(617)  377-3311 


MEASUREMENT  NEEDS  FOR  AURIC: 
ATMOSPHERIC  ULTRAVIOLET  RADIANCE  INTEGRATED  CODE 

OUTLINE 

Atmospheric  ultraviolet  radiation 
Comparisons  with  LOWTRAN  7 
AURIC  overview 


Measurements  needs 
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Atmospheric  Ullraviolel  Radiance  integrated  Code  (AURIC) 


MEASUREMENT  NEEDS  FOR  AURIC: 

ATMOSPHERIC  ULTRAVIOLET  RADIANCE  INTEGRATED  CODE 

MEASUREMENTS  NEEDS 

*  Atmospheric  density,  composition,  temperature,  etc. 

*  Cross-sections  for  absorption,  ionization,  etc. 

*  Reaction  rates  for  airglows 

*  Solar  flux  and  variability  (through  EUV) 

*  Transport,  scattering,  fluorescence 

*  Day/night;  Solar  cycle;  Geophysical  variability 

*  Auroral  regions,  particle  excitiation  processes 

*  Field,  laboratory,  associated  models 

*  Limb,  nadir,  temporal  variability 


MODELING  OF  ULTRAVIOLET  &  VISIBLE  TRANSMISSION  WITH  THE  UVTRAN  CODE 

E.  Patterson 

School  of  Geophysical  Sciences,  Georgia  Institute  of  Technology,  Atlanta,  GA  30332 

J.  Gillespie 

U.S.  Army  Atmospheric  Sciences  Laboratory,  WSMR,  NM  88002 

We  have  developed  a  model,  UVTRAN,  which  allows  the  calculation  of  transmittance  of  visible  and 
ultraviolet  wavelengths  as  well  as  the  calculation  of  lidar  returns  from  backscattcring  or  fluorescing  targets. 
The  transmission  model  is  designed  for  relatively  short  ranges  in  the  lower  troposphere  and  incorporates 
gaseous  scattering  and  absorption  as  well  as  aerosol  attenuation.  Of  special  note  is  the  treatment  of 
aerosol  attenuation.  The  aerosol  attenuation  is  parameterized  with  a  wavelength  dependence  that  depends 
0,1  the  visual  range.  The  usefulness  of  this  parameterization  relative  to  the  more  usual  parameterization  in 
terms  of  aerosol  microphysical  properties  whose  wavelength  dependence  is  independent  of  concentration 
will  be  discussed.  A  set  of  measurements  is  in  progress  that  is  designed  to  verify  the  transmittance 
models  will  also  be  discussed. 


MODELING  OF  ULTR/WIOLET  AND  VISIBLE 
TILVNSMISSION  USING  THE  UVTRAN  CODE 


by 


E.  M.  Patterson 

School  of  Earth  and  Atmospheric  Sciences  and  GTRI/EML 
Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332 

and 

J.  B.  Gillespie 

U.  S.  Army  Atmospheric  Sciences  Laboratory 
White  Sands  Missile  Range,  NM  88002 


WHY  UVTR/VN 


Simple,  user  friendly  model  for  transmission. 


Includes  major  molecular  absorbers. 


Direct  parameterization  for  wavelength  dependence. 


Modular,  suitable  for  inclusion  in  lidar  simulations  for 
both  DIAL  and  fluorescence  measurements. 


TRANSMISSION  MEASUREMENTS  TO  MODEL  VERIFICATION 


There  is  a  relatively  small  data  base  for  lower  tropospheric 
UV  transmission  along  horizontal  paths. 


Additional  measurements  under  a  variety  of  conditions  are 
needed  for  model  assessment  and  verification. 


ASL  has  a  program  for  model  verification: 

Urban  measurements  in  the  eastern  U.S. 

Co-operative  measurements  in  desert  areas  of  Middle  East. 
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IRANSMISSION  WORK 


AITENUATION  CALCULATIONS  FOR  STANDARD  CONDITIONS  AT  SEA 
LEVEL  FOR  AN  ASSUMED  VISIBILLIT  OF  10  KM  USING  THE 
UVTILAN  MODEL.  HIE  lOTAL  ATTENUATION  COEFFICIENT  IS 
(;IVEN  BY  THE  SOLID  LINE,  HIE  AEROSOL  AITENUATION  BY 
THE  DOTIED  LINE,  THE  MOLECULAR  SCATFERING  ATTENUATION 
BY  THE  DASHED  AND  DOITED  LINE,  THE  OZONE  AND  OXYGEN 
AITENUATIONS  BY  THE  LONG  AND  SHORT  DASHED  LINES, 
RESPECTIVELY, 


OBJECTIVES  OF  WORK 
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TRANSMISSOMETER  PATH  GEOMETRY 


ANCILLARY  MEASUREMENTS  NEEDED  FOR  MODEL  VERIFICATION 


Visual  Range 


Ozone 


Atmospheric  Pressure  and  Temperature 

Aerosol  Size-Number  Characterization 


SOITWARE  REQUIREMENTS 


CONTROL  i  ILTLR  WHEEL  POSITIONINO 

DATA  ACQUISITION 

INITIAL  DATA  PROCESSING 

CALCUALTION  OE  MEANS  AND  STD  DEVIATIONS 
EVALUATION  OE  DATA 

TRANSMISSION  CALCULATIONS 
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EXPERIMENTAL  CONSIDERATIONS 


CALIBItATION  NEEDS 

EFFECrrS  OF  UNEQUAL  PATH  LENGTHS 
BEAM  DIVERGENCE 

AMOUNT  OF  BEAM  SEEN  BY  RECEIVER 
ALIGNMENT  STABILITY 
SOURCE  STABILITY 

ATMOSPHERIC  REFRACTION  AND  TURBULENCE 


CALIBRATION  TECHNIQUES 


SIGNAL  IS  MEASURED  FOR  1TVO  PATH  LENGTHS  UNDER  HIGH 
VISIBILITY  CONDITIONS 


LONG  VVAVELENG1II  DATA  AT  TWO  DISTANCES  IS  USED  TO 
ESTIMATE  A  GEOMEIRICAL  CORRECTION  FACTOR 


SIMULTANEOUSLY  MEASURED  SCATTERING  CAN  BE  USED  TO 
CORRECT  LONG  WAVELENGTH  DATA 
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SUMMARY 


AUTOMATED  UV  TRANSMISSOMEl'ER  SYSTEM  HAS  BEEN  USED 
FOR  INirrAL  MEASUREMENTS  OF  ATMOSPHERIC  TRANSMITTANCE 
BEIAVEEN  200  AND  550  NM. 


FUI.UY  CAEIBR  Vn;!)  MEASUREMENTS  HAVE  NOT  YET  BEEN  MADE, 
BUT  IMIIAI.  DATA  INDICATE  THAT  SYSTEM  WILL  BE  SENSITIVE 
ro  VARIATIONS  IN  OZONE,  ATMOSPHERIC  DENSHT,  AND  AEROSOL 
(  HARACTFRISTK  S. 


SYSI  FM  IS  SI  I  I  ABI  i;  FOR  DESIRED  MODEL  V ERH  K'A HONS. 


SABLE:  A  SOUTH  ALANTIC  AEROSOL  BACKSCATTER  MEASUREMENT  PROGRAM 

S.B.  Alejandro,  G.G.  Koenig 
Geophysics  Laboratory/OPA,  Hanscom  AFB,  MA  01731 
J.M.  Vaughn,  P.H.  Davies 

Royal  Signal  &  Radar  Establishment,  St  Andrews  Road, 

Great  Malvern,  Worcs.,  WR143PS,  United  Kingdom 

The  plans  for  an  atmospheric  aerosol  field  measurement  experiment  are  described.  The  South 
Atlantic  Backscatter  Lidar  Experiment  (SABLE)  is  an  important  step  in  determining  the  feasibility  of 
obtaining  space-based  lidar  measurements  for  the  determination  of  atmospheric  winds,  and  the  spatial 
and  temporal  variation  of  aerosol  backscatter.  Two  field  campaigns  operating  from  Ascension  Island  (a 
British  colony  of  St.  Helena)  have  been  completed.  The  first  field  campaign  obtained  lidar  backscatter 
measurements  utilizing  an  airborne  platform  during  a  3-wcck  period  from  mid-October  through  mid- 
November  1988.  The  second  SABLE  program  was  expanded  to  include  measurements  utilizing  airborne 
particle-measurement  probes,  ground-based  lidars,  radiosondes,  and  supporting  meteorological  equipment. 
This  program  was  conducted  during  a  4-wcck  period  from  mid-June  Utrough  mid-July  1989.  This 
is  the  first  time  that  such  a  comprehensive  aerosol-measurement  program  has  been  undertaken  in  the 
Southern  Hemisphere.  SABLE-type  measurements  arc  not  only  pertinent  to  the  development  of  space- 
based  sensors,  but  also  have  direct  utility  in  the  study  of  aerosols  and  their  effect  on  climate. 


THE  RSRE  LASER  TRUE  AIRSPEED  SYSTEM  (LATAS) 

J.M.  Vaughn 

Royal  Signal  &  Radar  Establishment, 

St.  Andrews  Road,  Great  Malvern,  Worcs.,  WR143PS,  United  Kingdom 

I.A'I  AS  is  a  compact  robust  airlxtme  lidar  incoqiorating  a  10.6  /am  cw  CO2  waveguide  laser  of 
3  -4  Watt  output  and  currently  installed  in  a  Canberra  B-57  aircraft.  Doppler  shifted  radiation  scattered 
from  atmospheric  aerosols  is  detected  by  coherent  heterodyne  techniques.  Signal  processing  is  canned 
out  in  a  surface  accoustic  wave  spectrum  analyser,  A/D  converter  and  fast  integrator.  Integrated 
spectra  and  aircraft  parameters  are  recorded  on  a  high  density  MODAS  tape  recording  unit.  Detailed 
absolute  calibration  studies  have  been  carried  out,  together  with  extensive  algorithm  development  for 
data  extraction.  In  a  typical  flight  to  *15km  altitude  up  to  *10'*  individual  measurements  of  atmospheric 
back.rcaltcr  may  be  made.  Back.scattcr  sensitivity  extends  to  a  value  of  *8  x  lO^'^m  'sr'  as  currently 
operated.  Features  of  tlic  design,  calibration  and  perfomiance  will  be  discus.scd. 


MEASUREMENTS  OF  AEROSOL  AND  CLOUD  LIDAR 
BACKSCATTER  PROFILES  IN  THE  EQUATORIAL  SOUTH  ATLANTIC 
Lt.  Col.  G.G.  Koenig,  E.P.  Shcttlc*  and  S.B.  Alejandro 
Geophysics  LaboratoryADPA.  Hanscom  AFB,  MA  01731 
J.  M.  Vaughn 

Royal  Signal  &  Radar  Establishment, 

St.  Andrews  Road,  Great  Malvern,  Worcs.,  WR143PS,  United  Kingdom 

Aerosol,  cloud,  and  clean  air  lidar  backscatter  information  obtained  using  a  continuous  wave  coherent 
10.6  micrometer  CO2  laser  system  flown  on  a  B-57  Canberra  aircraft  will  be  presented.  These  data  were 
collected  during  the  South  Atlantic  Backsca'ter  Lidar  Experiment  (SABLE)  in  the  fall  of  1988.  SABLE 
is  a  cooperative  aerosol  lidar  backscatter  measurement  program  between  the  Royal  Signals  and  Radar 
Establishment  (RSRE)  of  the  United  Kingdom,  and  the  Geophysics  Laboratory  (AFSC)  of  the  United 
States.  The  program  was  organized  to  address  a  mutual  interest  in  establishing  a  database  of  atmospheric 
aerosol  backscatter  coefficients  in  and  around  a  remote  (Ascension)  island  location  in  the  Southern 
Hemisphere.  Of  special  interest  arc  the  lidar  signal  levels  associated  with  high  altitude  visible  and  sub- 
visible  cirrus.  Visible  and  sub-visible  cirrus  and  marine  stratocumulus  lidar  backscatter  coefficients  will 
be  presented.  The  observed  stratocumulus  lidar  backcattcr  coefficients  will  be  compared  with  model 
computed  backscatter  coefficient  information.  Finally,  the  backscatter  coefficients  associated  with  a 
documented  case  of  sub-visible  cirrus  will  be  described.  This  is  the  first  time  that  such  a  comprehensive 
aerosol  measurement  program  has  been  undertaken  in  a  remote  location  in  the  Southern  Hemisphere. 

*Now  at  NRL. 


7?. 7 


COMPARISONS  BETWEEN  SAGE  II  1.02  fim  EXTINCTION 
AND  SABLE  10.6  fim  BACKSCATTER  MEASUREMENTS 


E.P.  Shcttle*,  G.G.  Koenig,  and  S.B.  Alejandro 
GIVOPA,  Hanscom  AFB.  MA  01731 

J.  M.  Vaughn  and  D.W.  Brown 
Royal  Signal  &  Radar  Establishment,  Sl  Andrews  Road, 

Great  Malvern,  Worcs.,  WR143PS,  United  Kingdom 

The  South  Atlantic  Backscatter  Lidar  Experiment  (SABLE)  was  designed  to  provide  a  database  of 
atmospheric  aerosol  backscatter  coefficients  in  remote  sensing  of  the  South  Atlantic.  Such  data  would  be 
utilized  in  evaluating  the  feasibility  of  using  space-based  lidar  systems  to  obtain  global  measurements  of 
aerosols  and  winds.  Several  of  the  SABLE  flights  were  planned  to  give  temporal  and  spatial  coincidence 
with  SAGE  II  (Stratospheric  Aerosol  &  Gas  Experiment)  measurements  of  aerosol  extinction.  SAGE, 
which  is  a  satellite  based  instrument,  derives  the  aerosol  extinction  profiles  from  measurements  of  solar 
irradianccs  through  the  atmosphere  during  sunrise  and  sunset.  Comparisons  of  the  two  measurements 
provide  a  partial  validation  of  the  SABLE  measurements  and  a  validation  of  the  use  of  aerosol  models 
to  relate  the  SAGE  1.02  nm  extinction  to  the  10.6  /xm  backscatter  measurements.  These  aerosol  models 
can  then  be  employed  with  the  full  SAGE  aerosol  database  to  extend  the  spatial  and  temporal  coverage 
of  the  SABLE  measurements. 

*Now  at  NRL. 
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DEPARTMENT  OF  THE  AIR  FORCE 

GEOPHYSICS  LABORATORY  (AFSC) 

HANSCOM  AIR  FORCE  BASE.  MASSACHUSETTS  01731-5000 


REPLY  TO 

ATTN  OF;  QP  (F.X.  Kncizys,  617-377-3654/E.P.  Shettle,  617-377-3665) 

Annual  Review  Conference  on  Atmospheric  Transmission  Models 
TO  DISTRIBUTION  5  February  1990 


1.  The  extended  DoD  Plan  for  Atmospheric  Transmission  Research  and  Development 
tasks  the  Air  Force  to  conduct  an  annual  review  conference  to  provide  for  tri-service 
discussion  of  model  deficiencies  and  recommend  corrective  action. 

2.  The  13th  Annual  Review  Conference  will  be  held  at  GL,  Hanscom  Air  Force  Base, 
Bedford,  Massachusetts,  during  the  first  week  in  June  (5-7  June  1990).  The  objectives  of 
the  meeting  are  to  review  the  current  status  of  the  transmittance/ radiance  models  and 
ihe  need  and  plans  for  future  developments.  Areas  of  interest  include  molecular  and 
aerosols  effects  on  atmospheric  propagation,  propagation  measurements  using  lidar,  and 
turbulence  measurements  and  modeling.  We  are  interested  in  learning  about 
comparisons  between  models  and  experimental  data  and  in  hearing  your 
recommendations  regarding  how  we  can  overcome  any  identified  deficiencies. 

3.  This  letter  is  intended  to  solicit  contributions  to  this  conference.  If  you  would  like  to 
present  a  paper,  please  provide  an  unclassified  abstract  which  should  be  double-spaced 
and  no  more  than  12  lines.  Please  send  abstracts  to  G.P.  Anderson.  GL/OPE  (AFSO. 
Hamscom  AFB.  MA  01731-5000  bv  13  April  1990.  fTeiephone  617-377-23351. 

4.  Note  all  sessions  will  be  open,  and  it  is  anticipated  that  foreign  nationals  will  be 
permitted  to  attend  this  year.  Allow  sufficient  time  to  obtain  the  necessary  clearances 
for  your  papers. 

5.  If  you  plan  to  participate  in  the  conference,  but  not  make  a  presentation,  please 
notify  us  in  writing,  by  20  April  1990  (non-US  citizens  should  allow  at  least  6  weeks  for 
visit  approval).  Seating  accommodations  are  limited.  Send  your  letter  to  L.W.  Abreu, 
GL/OPE  (AFSC),  Hanscom  AFB,  MA  01731-5000,  (Telephone  617-377-2337).  If  you 
have  any  questions  contact  us  at  the  numbers  listed.  AUTOVON  prefix  for  Hanscom 
AFB  is  478. 


FRANCIS  X.  KNEIZYS 
Atmospheric  Effects  Branch 
Optical/Infrared  Technology  Div. 


I-RIC  P.  SHETl  LI-: 


Atmospheric  Optics  Branch 
Optical/Infrarcd  Technology  Div 
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DEPARTMENT  OF  THE  AIR  FORCE 
GEOPHYSICS  LABORATORY  (AFSC) 

HANSCOM  AIR  FORCE  BASE,  MASSACHUSETTS  01731-5000 


mplytoOP/F.  Kneizys/617-377-3654/E.  ShettlG/617-377-3665  1  May  90 

ATTN  Of 

Annual  Review  Conference  on  Atmospheric  Models 
su8jecT^^_g  9QJ  (AFGL/OP  Ltr,  5  Feb  90,  same  subject) 

DISTRIBUTION 

TO 

1.  On  the  basis  of  responses  to  our  February  5th  letter  announcing 
the  Annual  Review  Conference  on  Atmospheric  Transmission  Models,  we 
have  constructed  the  enclosed  tentative  agenda  for  a  two-day  meeting 
on  5-6  June  1990.  In  general,  we  have  reduced  some  initial  time  re¬ 
quests  slightly  to  allow  for  adequate  discussion.  Speakers  should 
plan  on  a  maximum  presentation  time  of  15  minutes.  Vugraphs  are 
preferable  as  visual  aids  in  the  talks.  The  meeting  will  be  un¬ 
classified  . 

2.  We  ask  that  all  conference  speakers  provide  us  with  hard  copies 
of  viewgraphs  used  in  your  presentation  at  the  meeting.  You  can  ei¬ 
ther  mail  them  to  F.X.  Kneizys  GL/OPE,  Hanscom  AFB,  MA  01731,  or 
bring  a  copy  to  the  meeting. 

3.  We  will  meet  on  5-6  June  (starting  at  0845  on  the  5th  and  0830 
on  the  6th  in  the  GL  Science  Center,  Bldg  1106,  Hanscom  AFB.  (See 
attached  map) .  Please  plan  to  arrive  by  0830  on  the  fifth  to  allow 
time  for  registration  on  the  first  day. 

4 .  We  do  not  intend  to  make  motel  reservations  for  any 
participants,  but  we  have  enclosed  a  list  of  motels  which  are  rea¬ 
sonably  close  to  Hanscom  AFB  (although  you  will  need  a  car  for 
transportation) .  On-base  government  quarters  are  limited,  so  if  you 
require  VOQ  accommodations,  please  make  your  reservations  early  by 
calling  617-377-2112.  The  AUTOVON  prefix  for  Hanscom  AFB  is  478. 

Francis  x.  kneizys  eric  p.  shettle 

Atmospheric  Effects  Branch  Atmospheric  Optics  Branch 

Optical/Inf rared  Technology  Div.  Optical / Inf rared  Technology  Div. 

3  Atchs 

1 .  Agenda 

2 .  Map 

3-  List  of  Motels 
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ANNUAL  REVIEW  CONFERENCE  ON  ATMOSPHERIC  TRANSMISSION  MODELS 

5-6  June  1990 

Geophysics  Laboratory  (AFSC) 

Hanscom  Air  Force  Base 
GL  Science  Center,  Building  1106 

PROGRAM 

Tuesday,  5  June  1990 
(0845  -  1200) 

WELCOME  -  Col.  Robert  J.  Hovde  (GL) 

KEYNOTE  -  Col.  T.S.  Cress,  USAF(OUSDA) 

SESSION  1  -  ATMOSPHERIC  PROPAGATION  MODELS 

(Co-Chair  Persons:  R.  Picard  (GL),  C.M.  Randall  (Aerospace),  R.G.  Isaacs  (AER)) 
“Status  of  the  LOWTRAN  and  MODTRAN  Models” 

L.W.  Abreu,  F.X.  Kneizys,  G.P.  Anderson,  E.P.  Shettlc,  J.H.  Chetwynd  (Geophysics  Laboratory) 

“Modeling  Solar  and  Infrared  Radiation  Fields  for  BTI/SWOE” 

J.R.  Hummel  (SPARTA,  Inc.) 

COFFEE  BREAK  (1020  -  1040) 

“PCTRAN  7  -  An  Implementation  of  the  GL’s  LOWTRAN  7  Model 
for  the  Personal  Computer” 

J.  Schroeder  (ONTAR  Corporation) 

“Spectral  Modeling  of  Off-axis  Leakage  Radiance  Using  the  MODTRAN  Code” 

N.  Grossbard  (Boston  College),  D.R.  Smith  (Geophysics  Laboratory) 

“Cloud  Opacity  Retrievals  Using  LOWTRAN  and  the  Stamnes  Scattering  Model” 

B.L.  Lindner,  R.G.  Isaacs  (Atmospheric  and  Environmental  Research,  Inc.) 

“Calculated  Cloud  Edge  SWIR  Radiance  Gradients  as  Viewed  by  Satellite” 

L.L.  Smith  (Grumman  Corporate  Research  Center) 


LUNCH  (1200  -  1330) 


Tuesday,  5  June  1990 
(1330  -  1700) 


(1330  -  1450) 

SESSION  2a  -  ATMOSPHERIC  PROPAGATION  MODELS  (CONTINUED) 

(Co-Chair  Persons:  R.  Picard  (GL),  C.M.  Randall  (Aerospace),  R.G.  Isaacs  (AER)) 


“The  Department  of  Energy  Initiative  on  Atmospheric  Radiation  Measurements 
(ARM):  A  Study  of  Radiative  Forcing  and  Feedbacks” 

R.G.  Ellingson  (University  of  Maryland),  G.M.  Stokes  (Pacific  Nortliwcst  Laboratories),  A.  Patrinos 
(US  Department  of  Energy) 


“The  HITRAN  Molecular  Database  in  1990” 

Lt.  S.  Shannon,  L.S.  Rothman  (Geophysics  Laboratory) 


“Line  Coupling  Calculations  for  infrared  Bands  of  Carbon  Dioxide  for  FASCOD3” 

M  L.  Hoke,  F.X,  Kneizys,  J.H.  Chelwynd  (Geophysics  Laboratory);  S.A.  Qougli  (Atmospheric  and 
Environmental  Research,  Inc.) 


“FASCOD3:  An  Update  (with  NLTE  Emphasis)” 

G.P.  Anderson,  F.X.  Kneizys,  E.P.  ShcttJc,  J.U.  Chetwynd,  L.W.  Abreu,  M  L.  Hoke  (Geophysics 
Laboratory);  S.A.  Gough,  R.D.  Worsham  (Atmospheric  and  Environmental  Research,  Inc.) 


COFFEE  BREAK  (1450  -  1510) 


Appendix 


Tuesday.  5  June  1990 
(1330  -  1700) 


(1510  -  1700) 

SESSION  2b  -  LINE-BY-LINE  APPLICATIONS 

(Co-ChaIr  Persons:  R.  Picard  (GL),  C.M.  Randall  (Aerospace),  R.G.  Isaacs  (AER)) 


“RAD:  A  Line-by-Line  Non-LTE  Radiative  Excitation  Model*' 

R.H.  Picard,  R.D.  Sharma,  J.R.  Winick  (Geophysics  Laboratory);  P.P.  Winterstcincr,  AJ.  Paboojian,  R.A- 
Joseph  (ARGON  Corp);  H.  Nebel  (Alfred  University) 


“Non-LTE  CO2  Vibrational-Temperature  Profiles  for  FASCOD3’’ 

P.P.  Winterstcincr,  A.J.  Paboojian  (ARGON  Corporation);  J.R.  Winick,  R.H.  Picard,  (Geophysics  Lab¬ 
oratory) 


“Non-LTE  CO  Infrared  Emission  In  the  Mesosphere  and  Lower 
Thermosphere  and  Its  Effect  on  Remote  Sensing” 

J.R.  Winick,  R.H.  Picard  (Geophysics  Laboratory);  P.P.  Winterstcincr  (ARGON  Corp.) 


“Path  Characterization  Algorithms  for  FASCODE” 

R.G.  Isaacs,  S.A.  Qough,  R.D.  Worsham,  J.-L.  Moncct,  B.L.  Lindner,  L.D.  Kaplan  (Atmospheric  and 
Environmental  Research,  Inc.) 


“Temperature  Retrievals  with  Simulated  SCRIBE  Radiances” 

J.-L.  Moncct.  S.A.  Clough,  R.D.  Worsham,  R.G.  Isaac.s,  L.D.  Kaplan  (Aimosphcric  and  Environmental 
Research,  Inc.) 


SOCIAL  HOUR  -  (1730  -  1900)  Officer’s  Club 


Appendix 


Wednesday,  6  June  1990 
(0830  -  1200) 


SESSION  3a  -  TURBULENCE 
(Chairperson:  J.H.  Brown  (GL)) 


“Validation  of  Random  Phase  Screens" 

A.E.  Naiman,  T.  Goldring  (W.J.  Schafer  Associates,  Inc.) 


“Recent  Advances  in  Modeling  of  Boundary  Layer  Refractivity  Turbulence” 

V.  Thcimiann,  A.  Kohnlc  (FGAN-Forschungsinslitul  fur  Optik) 


SESSION  3b  -  AEROSOLS 

(Co-Chair  Persons:  E.P.  Shettle  (NRL),  R.B.  Husar  (Wash.  U.)) 


“Visibility  Data  Filters  for  Europe” 

B.A.  Schichtcl,  R.B.  Husar  (Washington  University) 


COFFEE  BREAK  (0950  -  1010) 
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Appendix 


Wednesday.  6  June  1990 
(0830  -  1200) 


(1010  -  1200) 

SESSION  3b  -  AEROSOLS  (Continued) 

(Co-Chair  Persons:  E.P.  Shettle  (NRL),  R.B.  Husar  (Wash.  U.)) 


“Comparison  of  XSCALE  89  Vertical  Structure  Algorithm  with  Field  Measurements” 

R.P.  Fiegcl  (U.S.  Army  Atmospheric  Sciences  Laboratory) 


“Spatial  Non-Uniformity  of  Aerosols  at  the  15,500  ft  Level” 

L.A.  Mathews  (Naval  Weapons  Center).  P.L.  Walker  (Naval  Postgraduate  School) 


“A  Method  of  Estimating  Surface  Meteorological  Ranges  from 
Satellite  Aerosol  Optical  Depths” 

D.R.  Longtin,  J.R.  Hummel  (SPARTA.  Inc.);  E.P.  Shettle  (Geophysics  Laboratory) 


SESSION  4a  -  UV  MODELING 

(Co-Chair  Persons:  R.E.  Huffman  (GL),  J.B.  Gillespie  (ASL)) 


“A  Comparison  of  the  UVTRAN  and  LOWTRAN7  Aerosol  Models” 

S.  O’Brien  (Las  Cruces  Scientific  Coasulting) 


“High  Resolution  Solar  Spectrum  Between  2000  and  3100  Angstroms” 

I,. A.  Hall,  G.P.  Anderson  (Gcopliysics  Laboratory) 


LUNCH  (1200  -  1330) 


Appendix 


Wednesday,  6  June  1990 
(1330  -  1700) 


(1330  -  1510) 

SESSION  4a  -  UV  MODELING  (Continued) 
(Co-Chairpersons:  R.E.  Huffman  (GL),  J.B.  Gillespie,  (ASL)) 


“Predissociation  Line  Widths  of  the  Schumann-Runge  Absorption 
Bands  of  Some  isotopes  of  Oxygen” 

W.H.  Parkinson,  J.R.  Esmond,  D.E.  Freeman,  K.  Yoshino  (Harvand-Smithsonian  Center  for  Astrophysics); 
A.S-C.  Cheung,  S.L.  Chiu  (University  of  Hong  Kong) 


“AURIC:  Atmospheric  Ultraviolet  Radiance  Integrated  Code” 

R.L.  Huguenin,  M.A.  LeComptc  (Aerodyne  Research,  Inc.);  R.E.  Huffman  (Geophysics  Laboratory) 


"Measurement  Needs  for  AURIC:  Atmospheric  Ultraviolet 
Radiance  Integrated  Code” 

R.E.  Huffman  (Geophysics  Laboratory) 


“Modeling  of  Ultraviolet  and  Visible  Transmission  with  the  UVTRAN  Code” 

E.  Patterson  (Georgia  lastitute  of  Technology),  J.  Gillespie  (U.S.  Army  Atmospheric  Sciences  Laboratory) 


COFFEE  BREAK  (1510  -  1530) 
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Wednesday,  6  June  1990 
(1330  -  1700) 


(1530  -  1700) 

SESSION  4b  -  LIOAR 

(Co-Chairpersons:  Lt.  Col.  G.G.  Koenig,  S.B.  Alejandro  (GL)) 


“SABLE:  A  South  Atlantic  Aerosol  Backscatter  Measurement  Program” 

S.B.  Alejandro,  G.G.  Koenig  (Geophysics  Laboratory);  J.M.  Vaughn,  P.H.  Davies  (Royal  Signals  and 
Radar  Establishment) 


“The  RSRE  Laser  True  Airspeed  System  (LATAS)” 
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